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Preface 


The final report of work performed with support from NASA Grant NAG 
3-323 is presented in three parts. The first consists of three senior 
projects which examine: 

1) the sensitivity of the acousto-ultrasonic evaluation method to 
subtle forms of damage development in cyclically loaded composite 
material , 

2) the ability of the acousto-ultrasonic evaluation method for 
detecting initially existing imperfect material condition that promotes 
damage localization and final specimen failure, and 

3) the use of a well -characterized receiving sensor (an 
NBS/Proctor transducer) in performing acousto-ultrasonic evaluation of 
laminate composite materials. 

The second examines the nature of the wave propagation that occurs 
during acousto-ultrasonic evaluation of composite plates. It clearly 
demonstrates the part played by Lamb or plate wave modes and discusses 
how this may be utilized for characterizing such materials. 

Part III describes work done to replace the receiving transducer in 
the acousto-ultrasonic measurement procedure with calibrated noncontact 
“practical" optical detection system. 

A closure section presents a consolidated list of the conclusions 

of all these efforts. 
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Part I 


THREE SENIOR STUDENT RESEARCH PROJECTS 
INVOLVING THE CHARACTERIZATION OF COMPOSITE MATERIALS 
USING ACOUSTO-ULTRASONIC METHODS 


Seniors: Jonathan Buttram 

Paul Douglas 

Advisor: J.C. Duke, Jr. 

Senior: Thomas Lubnow 

Advisor: E.G. Henneke, II 


This work has been supported in part by NASA/LeRC through Grant 
NAG 3-323 


Materials Response Group 
Engineering Science and Mechanics 

Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 24061-4899 
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Preface 


The work herein is a result of work performed by three senior 

students in the Engineering Science and Mechanics Department in 

fulfillment of the requirement of the Senior Project and Report. The 

work is part of an ongoing overall effort to extend the understanding 
and increase the areas of application of Nondestructive Evaluation of 
Advanced Composite Materials by Means of Acousto-Ul trasonic Methods. 
Support in part for this work through the NASA/LeRC Grant NAG 3-323 
(technical monitor: Alex Vary) is gratefully acknowledged. The project 

advisors, J.C. Duke Jr. and E.G. Henneke, II, principal investigators of 
this grant, and Michael Kiernan, a graduate research assistant supported 
by this grant, have provided technical assistance. The financial 

support of the Engineering Science and Mechanics Department is 

gratefully acknowledged as well. 
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Assessment of Damage Development 
During Stiffness Degradation Plateaus 
in Cyclically Loaded Graphite Epoxy Laminates 


Paul B. Douglas 


Abstract 

The damage mechanisms that occur in cyclically loaded composite 
materials that experience a stiffness plateau are not known. This paper 
examines several nondestructive techniques which can be utilized to 
provide information on local and global material condition. It is found 
that the damage mechanisms affect sound attenuation in composite 
materials as shown by Stress Wave Analysis. 
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I. Introduction 


The increasing use of composite materials in critical structural 
components requires the designer to know specifics about their 
characteristics and behavior. Like composite materials' conventional 
counterparts, the "engineering metals", knowledge of the nature of the 
stress, strength/life, stiffness, and overall material condition is 
required to complete the design process and predict behavior. Unlike 
engineering metals, however, difficulties arise when defining composite 
characteristics and predicting behavior. First, composite properties 
are very direction-dependent; namely, individual lamina and laminates 
are anisotropic. Second, composites are the combination of two or more 
materials. As a result, lamina properties are a function of both the 

matrix and fiber properties. Lastly, as a consequence of both of the 
above, any behavior prediction process must consider the discontinuities 
and large gradients that arise at matrix-fiber interfaces and must 
include the favored failure modes and orientations. 

To complete composite material characterization, the previously 
mentioned inhomogeneities must be coupled with the element of time. In 
general, applied loads (stresses), ultimate strength or life, stiffness, 
and overall material condition are a function of both location in the 
material and time. Behavior prediction for composites, then, is very 
complicated and is a function of position and time. 

The goal of the design process and behavior prediction is to 
maximize component life and minimize cost and/or weight. Its ultimate 
goal is to accomplish this task and determine, or be able to determine 
through nondestructive methods, both the location and time component 
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performance will fall below design specifications -- when, where, and 
how failure will occur. 

Several schemes are available to predict failure in composite 
materials. The Tsai-Wu and Tsai-Hill equations are popular examples 
which rely on the introduction of ply stresses which are estimated by 
laminate analysis for a given load. In conjunction with laminate 
analysis, discount schemes are used to reduce stiffness in failed plies 
and permit calculation of new ply stresses as damage develops. This 
process of damage development and stiffness reduction is repeated until 
the failure criterion (equation) yields a load that falls below the 
design specification. 

Damage development in composite laminates subjected to cyclic 
loading is continuous. The degradation of stiffness, however, is not 
continuous and exhibits plateaus at various points. The damage 
mechanism(s) active during these plateau regions which leads to 
additional sitffness degradation are unknown. This makes it difficult 
to determine the condition of the material along the stiffness 
plateau. How close is the composite stiffness from falling off the 
plateau and subsequently failing? 

The objective of this project is to utilize laminates which exhibit 
long plateaus and apply nondestructive means of evaluation to identify 
the degradation mechani sm(s) . 
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II. Specimen Preparation 


This study was performed on T300/5208 graphite epoxy as supplied by 
NARMCO Materials, Inc. Specimens were fabricated locally with a flat 
face press using the cure time/temperature information found in Table 1 
[1]. The fiber orientation configuration was [0,±45) s . Specimens were 
cut using a circular saw with a diamond wheel blade. They measured 
1.021 by 10.0 inches and had a cured thickness of 0.032 inches. 
Specimen geometry is shown in Fig. 1. 

To improve gripping in the loading machine, the specimens were 
reinforced with Mi 1 ler-Stephenson Epoxy 907 on all surfaces for a length 
of 2.0 inches on each end. The edges were also reinforced with the 
epoxy except for the 1.0 inch gage section in the middle of the 
section. This was done in an attempt to concentrate damage development 
and enhance damage detection in the 1.0 inch gage length monitored by 
the extensometer. To accommodate the extensometer , aluminum tabs with 
V-shaped grooves were bonded to the specimen gage area with 3M 
Adhesive/Sealant Silicone (No. 03608). The reinforced specimen with 
attached extensometer tabs is shown in Fig. 2. 
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Table I. Cure Data 1 - for T300/5208 Graphite Epoxy Material 


Cure Temperature (°F) (°C) 275 (135) 

Cure Pressure ( ps i ) (MPa) 70 (0.5) 

Cure Time (min.) 105 

Post-cure Temperature (°F) (°C) 325 (163) 

Post-cure Pressure (psi) (MPa) 70 (0.5) 

Post-cure Time (min.) 120 

1 Supplied by manufacturer 
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90 


2. Specimen grip area and edge reinforcement, and extensometer V 
tab location diagram. 




III. Mechanical Testing 


Specimens were tested in tension-tension fatigue 

with a . / a = 0.1 and = 67.1 ksi [2j. This loading produces 

mm max max 3 r 

rapid initial degradation to the desired stiffness plateau. Fatigue 
tests were conducted using a MTS hyudraulic testing machine in the load 
control mode. Figure 3.0 is a block diagram of the testing machine and 
its monitoring equipment. Tests were run at 10 cycles/sec. Stiffness 
degradation was monitored by measuring the decrease in the secant 
modulus versus cycles. Strain measurements were made using an 
1.000 inch extensometer and were read directly from the MTS Digital 
Indicator in percent strain. 


IV. Nondestructive Testing 

The Stress Wave Factor (SWF)has been shown to be useful to detect 
the relative condition of a composite material. In tensile tests 
performed on composites, correlation between the point at which the 
specimen fails and the region having the lowest SWF value have been 
found [3]. Using this result, one can say that if failure occurs or is 
initiated at a point where maximum material degradation has taken place, 
then the SWF should be a means of monitoring relative material condition 
over time. The Stress Wave Factor is utilized to try to detect a change 
in the fatigue specimen condition which either shows the degradation 
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Fig. 3. Schematic of the mechanical testing and monitoring equipment. 






occurring during the stiffness plateau or marks the end of the plateau 
or both. 

During the cyclic tests. Stress Wave Factor values were recorded at 
points during the cyclic loading. Values were obtained at three (3) 
load levels -- minimum, mean, and maximum. The extensometer was removed 
from the specimen while these values were obtained. Two (2) one-half 
inch Panametric Transducers were secured to the specimen via a mounting 
block as shown in Fig. 4.0. Echo Laboratories' Ultragel II was the 
couplant used in the SWF tests. Equipment associated with the SWF is 
illustrated in Fig. 5.0. 

X-ray radiography is a method well suited for internal flaw 
detection in composite materials. Radiographic techniques can readily 
detect macroscopic matrix cracks and delamination in composite 
laminates. Specimens were prepared for each radiograph while still in 
the MTS loading machine. To enhance detection of flaws, a solution of 
60g zinc iodide, 10 ml water, 10ml isopropyl alcohol, and 10ml Kodak 
"Photo-Flo 600" was applied to the specimen before each radiograph. The 
solution was painted on the exposed edges of the laminate in the gage 
length while at mean load. After application, the specimen was cycled 
once manually using the set point load control, the zinc iodide solution 
was reapplied, and the cycle was repeated. The specimen was then 
unloaded and removed from the test machine. At this point, each 
specimen was washed with acetone and dried to ensure that none of the 
X-ray opaque zinc iodide solution remained on the surface. 

A Cabinet X-ray was used to make radiographs of the specimens. All 
radiographs wre made with an exposure of 25kV for three-tenths of one 
minute. Negatives were processed using a four-step developing 
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Transducer placement and spacing arrangement. 
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Schematic of the Stress Wave Factor testing system 
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Fig. 6. Normalized secant modulus versus loading cycles. (R-0.1). 


process: 1) 5 minutes in developer, 2) 30 second stop bath, 3) 3 
minutes in fixer, and 4) 10 minute water bath. All prints wre made on 
high contrast grade 5 Ilford black and white print paper. 


V. Results and Discussion 

The specified fatigue loading produced a rapid initial decrease in 
stiffness to the anticipated stiffness plateau. A graph of the 
normalized secant modulus versus cycles is plotted in Fig. 6.0. 

Specimen #1 was tested to 195,000 cycles at which time a decrease 
in stiffness was observed. This decrease, although seemingly small, 
occurred in a span of 10,000 cycles and was accelerating. The test was 
stopped and a visible edge delamination was found. It is believed that 
delamination was responsible for the drop off in stiffness. 

Stress Wave Factor values obtained using the procedure outlined in 
Section IV produced no correlation between the number of cycles and the 
degree of degradation of the specimen. Figure 7.0 shows this lack of 
correlation between stiffness and the SWF values. SWF values were taken 
at five (5) points during the 195,000 cycle test period of specimen 
#1. Obtaining values for these five (5) points involved stopping the 
loading machine after the required number of cycles. The transducers 
were attached as shown in Fig. 4.0 and the couplant was introduced to 
the transducer/specimen interface. The load was cycled several times to 
each test load using the set point, and SWF values were recorded. No 
general trend towards lower SWF values with respect to time at any of 
the test loads (min., mean, max.) could be found. In fact, 
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reproducibility of a SWF value at a given load for a given number of 
cycles was difficult. This difficulty in reproducing values is inherent 
in any type test involving the SWF. However, previous work has shown 
that by carefully controlling test variables, reproducibility of 10% is 
obtainable [3]. 

The test configuration described in Section IV did not permit good 
control of the test variables. Specific deficiencies were as follows: 

1* Vertical Test Position — allowed ultrasonic couplant to 
flow from between the transducers and test specimen under 
the force of gravity. The SWF value decreased with the 
length of time the transducers were left on the specimen. 

Transducer Placement/Orientation -- complicated the ability 
to "look" at the same location in the same way from point 
to point on the stiffness plateau. Either translation of 
the mounting block along the specimen or rotation of a 
transducer in the mounting block varied the SWF value. 

3. Elastic Connecting Bands — complicated the ability to 
reproduce transducer/specimen contact pressure. The SWF 
value decreased with increasing contact pressure. 

It was initially thought that obtaining SWF values at "load" would 
increase the detectability of material flaws. For instance, if flaws 
are some type of micro- or macroscopic crack, an applied load would tend 
to open the crack and increase sound impedance. This increase in 
detectability, however, did not account for, nor did it override, the 
loss of reproducibility imposed by the geometric constraints of the 
loading machine. 

Four (4) radiographs were made of specimen #1 during the 195,000 
cycle test. These are shown in Fig. 8.0. All except the last 
radiograph were taken when the material was still on the stiffness 
plateau. Figure 8.0d was taken after the decrease in the secant modulus 
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was observed. Also note that the V-tabs were removed from this specimen 
before the radiograph was made. 

First analysis of these photographs indicates that macroscopic 
damage is occurring while cycling along the stiffness plateau. It must 
be carefully noted, however, that damage (matrix cracks) is only visible 
because the zinc iodide solution has managed to work its way onto the 
cracks. Further, if the density of the photographed crack damage is 
related to the distance from the laminate edge (point of entry for the 
zinc iodide), a relation is seen. The material furthest from the source 
of the zinc iodide shows the lowest crack density. It should also be 
noted that the sequence of damage development shown in Fig. 8.0a-d would 
have most certainly produced a reduction in stiffness. With the above 
remarks noted, the following must be considered when analyzing zinc 
iodide enhanced X-ray radiographs of cyclically loaded composite 
materials: 

1. A radiograph does not show damage; it indicates where the 
enhancing dye is present. 

2. A sequence of radiographs does not show crack development 
but rather the extent of dye penetration. 

3. As a result of 1. and 2., and by observing the sequence of 
photrographs in Fig. 8.0, it can be proposed that the 
degree of macroscopic cracking between two (2) points on 
the stiffness plateau is very nearly the same, but might 
appear quite different with comparison of radiographs taken 
at each point. 

4. The cyclic nature of the load may act as a pumping 
mechanism to get the dye (zinc iodide) into the developed 
cracks. 

A somewhat different approach was used for the analysis of specimen 
#2. Seeing that the Stress Wave Factor (SWF) data for specimen #1 
yielded no correlation between damage development and position along the 
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stiffness plateau (cycles), it was decided to modify the signal 
interpretation. Increased sensitivity and reproducibility were 
accomplished by removing the testing deficiencies imposed by the loading 
machine. Namely, SWF data was obtained with the specimen removed from 
loading machine grips and positioned on a horizontal surface. The 
transducers were mounted as shown in Fig. 9.0 and their initial position 
was marked on the specimen with a white marker. Signal interpretation 
was changed from recording the peak freguency count to performing a Fast 
Fourier Transform of the Stress Wave Signal and looking at the zeroth 
moment" of this transform (5). The zeroth moment is the square root of 
the area under the frequency versus amplitude curve and gives a measure 
of the energy reaching the receiving transducer. Plots of the peak 
frequency count (amplitude versus time) and the transformed Stress Wave 
Signal (frequency versus amplitude) are shown in Fig. 10.0 and Fig. 11.0 
respectively. 

The zeroth moment calculation was performed at three (3) points 
along the stiffness plateau curve of specimen #2 plotted in Fig. 12.0. 
Preliminary results shown in Fig. 12.0 correlate a decrease in the 
moment calculation with increasing cycles. This is a very significant 
result because it indicates a continuous degradation of material 
condition not associated with a change in stiffness. As the composite 
was cycled and damage progressed, there was an increase in signal 
attenuation. 

Two (2) radiographs of specimen #2 were made. They are shown in 
Fig. 13.0. The comments made concerning radiographic evaluation of 
fatigued composites from the discussion of specimen #1 apply. 
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Testing configuration for stress wave analysis with the speci 
testing machine. 
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Fig. 10. Time domain received signal as gated for analysis. 
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Fig. 11. Fast Fourier Transform (FFT) analysis of the time domain 
signal of Fig. 10. 
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Fig. 13. Penetrant enhanced X-ray radiographs of specimen #2 after a) 70,000, b) 170,000 
loading cycles. 


VI. Conclusions 


The following conclusions are made: 


1. Zeroth moment calculations indicate an increase in sound 
attenuation with increasing cycles in a [0,±45 ] s composite 
laminate, with no corresponding increase in observable 
damage. 

2. Stress Wave Frequency Counts do not correlate with the 
continuous degradation occurring in a [0,±45l s composite 
laminate. 

3. X-ray radiography yields misleading visual information 
about the degree of material degradation occurring in a 
[ o ,±45 ] s composite laminate under cyclic loading. 
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A Study of the Final Failure Event 
for Uniaxial Tensile Loading 
of Graphite Epoxy Laminates 


Jonathan D. Buttram 


Abstract 

Identifying areas of low strength in a composite material is a 
problem of which little is known. This experiment approaches this 

problem by identifying damage characteristics that are related to areas 
of low strength through the use of stress wave factor analysis, edge 

replication, and X-ray radiography. 

The data obtained indicated a relationship between areas of low 
structural strength and those which produced low stress wave factor 
(SWF) values. Using this relationship, which provided the means to 
predict failure areas, and the ability to visually monitor these areas 
through the use of replicas and radiographs during increasing load 

increments, data was obtained which suggests that damage is a local 
phenomena that eventually governs failure. Visual damage 

characteristics found in the edge replicas and radiographs were found 
not to be as significant as stress wave factor data, which, if used 

properly, can be a characteristic of otherwise undetected damage. 
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I. Introduction 


The identification of potentially weak areas in a composite 
material has become a major problem in the forecasting of durability. 
Although there are many nondestructive testing techniques that can be 
used on composites, the interpretation of the results from these tests 
can be difficult and unsure. If certain distinguishing damage 
characteristics were to be found which correlated directly to areas of 
structural failure, then these could be used to assist in the 
nondestructive evaluation process. This project addresses this problem 
by initiating damage in a quasi-isotropic material by applying uniaxial 
tension to sample specimens. Using nondestructive testing techniques 
such as X-ray radiography, acoustic emission monitoring, stress wave 
factor analysis, and edge replication, a direct relationship between 
some distinct damage feature and the probable failure area is sought. 
To understand this more clearly, certain assumptions must be discussed 
and related to the procedures and final outcome. 

First one must ask the question, "Why does a failure occur at a 
particular location?" Theoretically if a material were perfect with no 
flaws, then a uniform specimen would just disintegrate under uniaxial 
tension. Since every point on the specimen has the same strength 
properties with no defects, then failure should simultaneously occur at 
every point along the specimen. This of course does not happen in 
reality because of imperfections in the material. In this experiment 
the assumption is made that as the specimen becomes stressed under load, 
the point with the most significant damage becomes increasingly more 
damaged. Immediately before failure the damage here is more advanced 
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than at other points throughout the sample. The density of this damage 
is also assumed to decrease as one moves farther from the failure 
area. Therefore if long specimens are used and the approximate location 
of failure can be determined, then it is assumed that another area that 
is much farther away will be less damaged. 

The detection of the failure area presents another problem. 
Because of the dynamic nature of failure under high stresses, unwanted 
destruction can be imposed on the material if actual rupture occurs 
during testing. Because of this, a stress wave factor analysis was used 
to predict areas of failure before failure occurred. Stress wave factor 
(SWF) is a relatively new nondestructive testing technique which 
utilizes high frequency stress waves to predict relative material 
strengths [2,3], This method assumes that as a stress wave propagates 
through a material, imperfections such as voids, cracks, defective 
bonds, broken fibers, etc., will attenuate the signal. During the 
testing process, high frequency pulses are transmitted by one transducer 
and received by another that is at a fixed distance away. A scanning 
process is set-up so that an entire specimen can be examined, a segment 
at a time. Accordingly, if a segment of a material exhibits a signal of 
particularly low amplitude, then one can assume the presence of 
imperfections or areas of low strength. This method is vital in this 
experiment since it allows the prediction of low strength areas that can 
be scrutinized during the testing and analysis process. 

II. Objective 

The purpose of this experiment is to find distinguishing damage 
characteristics that are directly related to predicted areas of low 
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structural strength through the use of acoustic emission, stress wave 
factor, edge replication, and X-ray radiography. 

Approach 

The basic experimental approach was to use the stress wave factor 
technique to first predict an area of low structural strength. Then 
samples were loaded to various levels and radiographs plus edge replicas 
were used to observe internal damage. The objective was to unload the 
specimen immediately before failure was expected, paying particular 
attention to the area of low strength predicted by the SWF analysis. 
Acoustic emission monitoring was also used so that failure could be 
anticipated from emission occurring during loading. 

III. Procedures 

A. Sample Preparation 

Six samples were cut from a [ + 45,-45,90,0) s , a graphite epoxy 
laminated sheet. These samples varied from 12 to 16 inches long and 
were H inch wide. Because of their length, some of the longer samples 
were cut so that they would fit the testing machine. All specimens were 
cut using a diamond impregnated cut off wheel employing continuous 

coolant flow. 

B. Testing Platform Preparation 

To achieve consistent results from the stress wave factor 
analysis, a position platform was designed (Fig. 1). The platform held 
the sample securely and provided a means by which to position the 
transducers accurately in the center of the specimen. Quarter inch marks 
were placed on the transducer rest so that a precise location of the 
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specimen could be recorded and reproduced. The transducers were 
maintained 1 inch apart by a simple fixture also shown in Fig. 1. The 
previously mentioned tools were used in the stress wave factor analysis 
in order to improve the accuracy of the results. 

C. Stress Wave Factor 

An understanding of the basic equipment set-up and its 
function must be obtained before the procedures for the SWF can be 
appreciated. The first step in the process is the production of 
ultrasonic pulses. These pulses are produced by an ultrasonic pulser 
which can be adjusted depending on the desired energy and repetition 
rate desired (Fig. 2). The pulse then travels to the transmitting 

transducer and propagates through the specimen to the receiving 
transducer. The signal is then amplified, filtered, and gated. This 
gated signal is sent to a counter which counts cycles over an adjustable 
threshold (voltage level) in a desired time period. The number of 
cycles counted was recorded according to specimen location. In order to 
smooth out any abrupt changes in data values, a computer program was 
written. This program took the data values and averaged each with its 
two neighbors. For example: 

i (Xp + *3 + X 4^ 

X 3 = O 


or 


(X„-l + X n + x n + l> 


3.0 
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Fig. 2. Schematic of the stress wave factor testing set-up. 


where X n represents the actual data and X' n represents the new 
calculated data. This method only made data more plottable but did not 
change data trends which are essential in this experiment. For more 
detailed information pertaining to topics concerning SWF, refer to the 
appendix. 

The following steps were preformed in the analysis. 

Step 1: The specimen was placed in the positioning platform and 
adjustments performed so that the center of the transducers were 
positioned directly above the center (of width dimension) of the 
specimen. A small overhang of the transducers was experienced in this 
case because of the dimensions of the transducer and specimen. 

Step 2: The transmitting and receiving transducers were fixed at 1 
inch apart. A starting point was designated and a large drop of 
couplant was applied, and another, 1 inch along the sample. The 
couplant was then spread into a thin, even film using a razor blade. 

Step 3: The transducers were set in the couplant and a 50 gram 
mass was placed on each transducer. The transducers were allowed to set 
for 25 seconds so that an even and uniform layer of couplant was 
achieved. 

Step A: After the 25 seconds had expired, the counter was reset 
and signal preamplifier turned on. Twenty more seconds were allowed to 
pass and a reading from the counter was then recorded. 

Step 5: The transducers were removed from the specimen and all 
couplant cleaned from the transducer and specimen surfaces. 

Step 6: The tandem transducer assemply was translated H inch 
along the specimen to the next measurement location and steps 2-5 
repeated until the entire testing area was examined. (Two inches on the 
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ends of each sample were not tested because they were located in the 
testing machines grips.) 

Step 7: The data was then entered into a computer program which 

calculated and plotted the results. 

0. Sample Test 

A preliminary test was run on one of the samples in order to 
determine a representative ultimate strength. 

E. Edge Replication Technique 

Edge replication is a simple and quick method to examine 
damage on an exterior surface. It provides a "f ingerprint" of the 
desired surface showing damage such as cracks, delaminations, voids, 
etc. A special tape is used which softens with the application of 
acetone. If pressed against a surface during this soft condition, the 
tape will be impressed with features from that surface. The steps in 
the procedure are as follows. 

Step 1: The edges of the sample to be tested were cleaned. (These 

surfaces should be smooth.) 

Step 2: A thin piece of replica tape was cut from the roll and 

securely fastened to the edge of the specimen by taping it at the ends. 

Step 3: Acetone was then applied to the crevice between the 

replica tape and specimen edge by means of a syringe. 

Step 4: An object with a flat, smooth surface was rubbed over the 

back of the tape, pressing it against the sample surface. Approximately 
45 seconds was allowed to pass before the tape was removed. This 

provided time for the acetone to dry and the tape to harden. 

Step 5: The samples were marked where replicas were taken. The 

edge replicas were placed on a glass slide and viewed under 
magnification. 
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F. Second Stress Wave Factor 

Another stress wave factor analysis was performed after load 
of the specimens. (Refer to part C for the precise procedure.) These 
results were later compared to the SWF results obtained prior to loading 
so that a more accurate prediction of the low strength area could be 
made. It was believed that the 1000 lb. loading might produce damage in 
the predicted weak spots and enhance the SWF data. 

G. Tests to Failure 

Step 1: Each specimen was positioned in the Instron constant 
crosshead rate testing machine grips so that uniaxial tension was 
applied with no shear forces or bending moments present. A transducer 
was clamped to the specimen which monitored the acoustical signal. This 
signal was directed into a preamplifier and then to a high pass 
filter. After being filtered, the signal was further amplified and sent 
to an oscilloscope that provided a real-time display. This set-up (as 
shown in Fig. 3) was used during the monitoring of the acoustic signals 
emitted during the damage process. 

Step 2: The sample was then loaded to 1200 lbs. and zinc iodide 
(Znl 2 was generously applied to the samples' edges. The Znl 2 was 
allowed 3-4 minutes to penetrate all the cracks that were present. The 
Znl 2 , which is radiographically opaque, was applied so that any 
delaminations or cracks would be enhanced on the X-ray radiographs. 
Edge replicas were then taken and the specimen was unloaded. X-ray 
radiographs were made. 

Step 3: The sample was reloaded to 1250 lbs again following the 
procedures outlined in step 2. 
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Step 4: 


The sample was loaded to 1300 lbs. The procedures 

explained in step 2 were repeated. 

Step 5: The specimen was reloaded and the load increased again in 

50 lb. increments while the procedures explained in step 2 were repeated 
until either the breaking of zero degree fibers was heard, the next 
increment was reached, or failure occurred. Because of the high 
stresses incurred at this stage, the specimens were unloaded to 
approximately 900 lbs. for the application of Znl 2 ^ replication 

procedure. This was done so that further damge would not occur while 

the load was maintained. 

H. X-ray Radiography Technique 

X-ray radiographs were taken so that any internal flaws could 
be detected. A cabinet type X-ray machine was used in this 
experiment. A tube voltage of 25kV was used with an exposure time of .9 
minutes. 


IV. Results 

Six samples were tested in this study. Three of the specimens 

broke in or near the testing machine's grips. Knowing that the grips 
can cause undesired stresses, initial nondestructive examination was not 
performed on these areas. Therefore these specimens were not given 
further consideration because of the lack of SWF, X-ray and replica 
data. The location of failure and four lowest SWF measurements are 

tabulated below. 
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Position Predicted by SWF Position Predicted by SWF 


Sample 

Number 

Position 

of 

Failure 

(inches) 

Before Failure (inches) 
12 3 4 

After Loading (inches) 
12 3 4 

1 

GRIPS 

— 

— 

— 

— 

— 

-- 

— 


2 

11.25 

12.0 

12.25 

12.5 

11.75 

8.5 

6.5 

8.75 

12.35 

3 

4.0 

5.75 

6.0 

5.5 

3.0 

4.25 

4.5 

3.75 

3.0 

4 

6.0 

0.5 

0.25 

6.5 

6.75 

7.25 

7.0 

.25 

1.0 

5 

GRIPS 

— 

— 

— 

— 

— 

— 

— 

— 

6 

GRIPS 

__ 









V. Data Analysis 

Stress Wave Factor 

When examining the data from the stress wave factor tests (Figs. 
4,5,6), a number of observations can be made. The first of these deals 
with the correlation between the predicted site of failure and the 
actual rupture location. Remembering that the low strength area should, 
by assumption, correspond to the lowest valued data points, a comparison 
between the failure site and the obtained data can be made. 

Sample 2 broke at position 45 which was located 11.25 inches from 
the end of the sample. As shown in Fig. 4, the SWF before the 1000 lb. 
loading predicted an area of low strength at 12.0 inches. If the 4 
lowest points are considered, then a more accurate 11.75 inches can be 
obtained. The SWF fater the specimen had been subjected to the 1000 lb. 
loading predicted (lowest point) that failure should have occurred at 
8.5 inches (Fig. 5). This location corresponds to the lowest point, 
with 12.25 inches being the fourth lowest and most accurate. This seems 
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Normalized SWF versus specimen position for specimen #2 before tensile loading. 
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Normalized SWF versus specimen position for specimen #2 after loading but before 
fai 1 ure. 


to suggest that failure can occur at low data regions that are not 
always the lowest in value. This topic will be discussed further, later 
in this paper. 

Examination of sample 3 showed that failure occurred 4 inches from 
the end of the sample. The stress wave factor analysis before and after 
the 1000 lb. loading indicated areas of low structural strength at 5.75 
and 4.25 inches from the end, respectfully. These again refer to the 
lowest values obtained from the SWF data. Table 2 illustrates that more 
accuracy can be obtained if more low points are considered. Plots of 
the SWF for sample 3 are found at Figs. 6 and 7. 

The last sample, 4, failed at exactly 6 inches from the end of the 
sample. Figure 4 illustrates that the data obtained before the 1000 lb. 
loading predicted a weak area at 6.5 inches (Fig. 8). The data obtained 
after the 1000 lb. loading (Fig. 9) indicates a low area at 7.0 inches 
from the end. 

The results above indicate remarkable accuracy if certain 
limitations are realized. Observations revealed that the procedures 
used entailed enough experimental error that the lowest value obtained 
does not necessarily correspond to the position of failure. A more 
accurate method for analyzing this data would be to look at the relative 
low spots on a normalized plot and be cautious of these areas. For 
clarity, consider the chart on the following page. 
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LOW POINTS 


ERROR (inches) 


CONSIDERED 

BEFORE 

LOADING 

AFTER 

LOADING 

1 

2.67 

1.5 

2 

2.67 

1.41 

3 

.92 

1.25 

4 

.67 

.75 


Table 2. Averaged Error Values Between Predicted and Actual Failure 
Positions. 
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Normalized SWF versus specimen position for specimen #3 before loading. 
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Normalized SWF versus specimen position for specimen #3 after loading. 
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Normalized SWF versus specimen position for spec 


X 
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POSITION < inches) 


This chart illustrates the averaged distance that the predicted 
position differed from the actual position of failure if the indicated 
number of low spots were considered. For instance, if the SWF data 
obtained after the 1000 lb. loading was used and if the positions of the 
four points of lowest value were considered, the failure occurred 
approximately .75 inches from one of these points. As the number of low 
points considered decreases, so does the location accuracy. 

Examining the SWF vs. POSITION plots given in this section, a 
tendency for low valued data points to bunch up in certain areas is 
noticed. Therefore if 4 low points are considered, and 3 of those 
points are located in one general area, then only two areas would have 
to be concentrated on for failure. The error related to this case would 
then only correspond to two locations instead of four, making the 
overall total area of predicted failure less. 

Another interesting observation that is found in the previous chart 
is that the predictions made after the 1000 lb. loading appear better 
than those taken before the loading if only the lowest two points are 
considered. As more low points are considered, the two accuracy values 
become much closer in value. This seems to illustrate that the SWF data 
obtained before loading had low areas that became even lower relative to 
their neighbors after loading. In other words, areas which were weak 
but not the weakest before loading were transformed into the weakest 
positions during the initiation of damage. Since the first SWF analysis 
predicted areas of low structural strength, then these areas should be 
the first to develop damage during the 1000 lb. loading, causing the 
second SWF to become more pronounced. Figures 10, 11, and 12 are 
superimposed graphs of the before and after loading SWF plots. These 
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NORMALIZED 
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versus specimen position for specimen 
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POSITION C inches > 


graphs show that in most cases the low spots remain low or become 
lower. There are cases where previously low spots become high, but 
these are not very common. Considering the numerous experimental 
variables that could affect the outcome (see error section), the b efore 
and after loading SWF plots are basically similar. The accuracy of the 
SWF data was estimated to be 5% based upon previous experimental tests. 

Edge Replicas 

Figures 13, 14, and 15 display prints of replicas taken of the 
failure areas prior to actual failure. They all show damage that 
eventually led to the failure. 

Sample 2 shows a severe longitudinal crack which reaches the edge 
of the sample at the location of failure. The crack is believed to have 
started at the failure area and then propagated down through the 
sample. This damage appears to be a combination of transverse cracks 
and delaminations. The crack uses the 90 deg ply to propagate for most 
of its length, except near the failure area where it breaks through the 
+45, -45 plies. 

Sample 3, Fig. 14, again shows delamination associated with its 
failure region. The print illustrates a crack that is relatively smooth 
with no abrupt changes in the direction until the failure area is 
reached. Here small, transverse cracks are present. These are believed 
to have grown, eventually passing through the material to the edge, 
causing failure. As in sample 2, damage is mostly restricted to the 
90 deg plies until the failure area is reached, where the +45, -45 plies 
are damaged. 
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Failure 

location 



Failure 

location 


a) 


b) 


ig. 13. Edge replicas of specimen #2 after loading close to the 
failure load, a) 20X b) 40X magnification. 
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Failure 

location 


a) 




Failure 

location 


o) 


Fig. 14. Edge replicas of specimen #3 after loading close to the 
failure load, a) 20X b) 40X magnification. 
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Sample 4, Fig. 15, indicates more localized type damage. Very 
little delamination can be found in this sample. Cracks are found to 
have passed through the 90 deg ply and started into the +45, -45 
plies. The damage is found on both sides of the 0 deg ply, which does 
not occur in the previous two samples. 

When comparing all three samples, it is observed that damage does 
not occur in the 0 deg ply. Therefore it can be said that the outer 
plies play a major role in the failure process. Failure of the outer 
plies causes an increase in the amount of stress applied to the 0 deg 
plies in a particular area. This increased stress leads to a failure at 
this area. Another question to be raised deals with the sequence of 
events during the failure process. In studying the replicas, one must 
conclude that the damage starts in the 45 deg plies and works its way to 
the surface by means of delamination and transverse cracking. 

X-ray Radiographs 

The radiographs shown in Fig. 16 are of sample 2 and confirm the 
delamination shown by the replicas. The radiographs are for increasing 
load increments. 

The first radiograph (low load) indicates a higher density of 90 
deg matrix splits at or near the failure area. Although the entire 
sample is not shown here, the damage shown is by far the most visible. 
This initial damage must play a role in the development of damage 
located on the opposite side of the sample, as shown in the next 
radiograph (medium load). This radiograph indicates a delamination 
which starts at the edge plus a large number of 90 deg matrix cracks. 
These cracks are observed throughout the material and not just at the 
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location 



a) 



Fig. 15. Edge replicas of specimen 44 after loading close to the 
failure load, a) 20X b) 40X magni fication. 
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a) b) 


Fig. 16. Penetrant enhanced X-ray radiographs at an a)low 
c) high level of tensile load. 


: : fif 



b) medium 
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eventual rupture location. The initial damage located on the opposite 
side of the sample had not changed very much at this point. The last 
radiograph shows how the area of delamination grew even more with 
increased load (approximately 1.5 times as big). This radiogj'raph was 
taken very close to actual failure and corresponds to the edge replica 
for this sample. 


VI. Error Analysis 

The purpose of this section is to briefly examine the methods of 
research and to reach a determination of the accuracy and reliability of 
the resulting data. 

Stress wave factor methods took many hours of trial and error 
sample runs to develop. The procedure which was developed was strictly 
followed with regard to variables such as transducer positioning, 
couplant application, instrument setting variations, transducer weights, 
surface preparation, etc., being carefully controlled. For example, 
every event in the testing procedure that could be time dependent was 
carefully timed. A special testing platform was designed to control 
transducer positioning. To control couplant film variations, practice 
runs were performed until a smooth consistent film could be applied. To 
assure adequate preamplifier battery power throughout the testing period 
for each sample (approximately 2 hrs.), the preamplifier was turned on 
only 20 seconds for each measurement. Of this 20 seconds, the first 10 
provided time for any transient effects to dissipate and the last 10 
seconds were used for testing. To control unclean surface problems the 
transducers and sample surface were cleaned after every measurement. To 
assure good transducer to sample contact, special weights were used. 
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Other controls, such as using the same transducers to transmit and 
receive each time, insured reasonably reproducible results. The error 
in this process was estimated to be approximately 5% and markings 
indicating this have been placed on all graphs in this report. 

The accuracy of the edge replicas should be high. Errors dealing 
with the replicas come basically with their examination. One must not 
mistake dirt, dust, fingerprints, etc., for flaws in the material. Good 
replicas are a product of practice and patience from which comes 
accurate prints. 

The X-ray radiographs can sometimes be misleading, so caution 
should be taken when examining them. Zinc iodide was used in the 
radiographic procedures so that internal flaws would appear on the 
film. It was applied while the samples were under load so that it could 
seep into small cracks. The problem here is that the zinc iodide does 
not penetrate completely into the material. Because of this the center 
of the sample may appear undamaged; this seems to be the case in this 
project. No damage appears to have occurred in the center of the 
samples on the radiographs. Whether or not this is so cannot be 
determined with the techniques used here. Therefore all conclusions 
drawn from the radiographs are most appropriate for the outer edges and 
not for the center of the sample. 

The next question that must be raised deals with the reliability of 
data and resulting conclusions. Since only three samples were 
considered for analysis, a reliability problem exists. With such a 
small sample number, the overall reliability cannot be determined. 
Although the results from the project may possibly be accurate, their 
reliability cannot be determined without further testing. Three tested 
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samples do not lay a foundation for a sufficient statistical basis. 
Because of this fact, the confidence level of these results is low. 
Therefore the conclusions from this experiment should not be used for 
the basis of any engineering decisions, but rather be used to motivate 
further investigation. 


VII. Conclusion 

It can be concluded from the results produced by this experiment 
that there are damage characteristics that are directly related to areas 
of failure in a graphite epoxy composite laminate under uniaxial 
tension. The edge replication and X-ray radiographic techniques were 
good detectors of damage after it had been established in the 
material. The material used in this experiment showed delamination and 
transverse matrix cracks that started and grew in the failure area. In 
the samples tested, these distinguishing features seemed to be good 
characteristics of the failure area since their severity increased in 
these areas. Edge replication and X-ray radiography were very useful in 
locating a failure location after gross amounts of structural damage had 
already been initiated. Visual techniques require an extensive amount 
of damage just to be detected. What is to be done if the damage is not 
visual and a failure area is to be established? To answer this question 

the SWF results must be interpreted. 

The results from the stress wave factor analysis point to the fact 
that damage can be detected. The results showed that predictions of 
locations of low strength areas in a material can be made with an 
accuracy of less than ±3/4 inch within the limitations discussed in the 
results section of this paper. Since structural damage in a material 
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affects its strength, and low SWF values seem to be related to low 
strength areas, structural damage should be related to low SWF values. 
Based upon this reasoning, one can say that low SWF values are a damage 
characteristic. It should also be pointed out that imperfections such 
as those encountered during the production process of the material can 
also cause low SWF values. Because of this, a low SWF values does not 
necessarily mean that damage is present. 

Another interesting conclusion appeared while examining the 
results. Because of the fact that visual damage was observed in the 
same general location of low SWF values and this damage was the first 
significant damage to appear on the replicas and radiographs, it can be 
said that initial damage is a local event. In other words, the damage 
does not initiate as a global event, occurring everywhere on the sample 
and then localizing immediately before failure occurs. This experiment 
suggests that initial damage is a local event and although damage occurs 
throughout the sample during further loadings, the failure will occur at 
this localized area. 
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Appendix 

Instrument Specifications 


A. Stress Wave Factor Set-Up: 
Transducers 

Preamp 

Pulser & Gate 
Counter 


Panametric Inc. 

Width 11/16 in. 

Freq. approx. 5 MHz 
Serial no. 43026 
43022 

Panametrics Ultrasonic Preamp 
40-60 db. 

Panametrics Model 5052ua 
Ultrasonic Analyzer 

Hewlett Packard 5326 Timer- 
Counter-DVM 


B. Instrument Settings for SWF: 


Panametrics Ultrasonic Analyzer 
Rep. Rate 
Damping 
Energy 
Attenuation 


200 

200 

2 

6 db. 


Width 

Length 

Thickness 


Sample Specifications 

8/16 in. 
12-16 in. 
.041 in. 
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The Use of an Improved Transducer 
For the Study of Wave Propagation in 
Composite Plates 


Thomas S. Lubnow, Jr. 


Abstract 

The performance of the Proctor Conical Transducer developed at the 
National Bureau of Standards as applied to composite materials has been 
evaluated and deemed satisfactory with the exception of the method of 
electrical grounding employed. Incorporation of the conical transducer 
into the Stress Wave Factor scheme yielded evidence of the possible 
plate wave mode of propagation in thin composite plates. 
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I. Introduction 


Fiber-reinforced composite materials are gaining widespread use in 
industrial and military settings. Unfortunately, the development in 
understanding of mechanical behavior and of failure theories has not 
kept pace with the' rapid increase in applications. This is in part 
because delaminations, voids, matrix cracks, and broken fibers 
contribute to a very complex damage state in composites. An increased 
dependence on nondestructive testing and the development of new NDT 
techniques is a result of the need to better understand the performance 
of these materials. One such technique is termed acousto-ul trasonic 

testing or acousto-ultrasonics (AU). 

A relatively recent development in the field of nondestructive 
evaluation of materials, acousto-ultrasonics combines the advantageous 
aspects of acoustic emission and ultrasonic testing techniques. With 
this technique, the acoustic pulse is generated outside the specimen, 
sent through the material, and received by another transducer a short 
distance away (Fig. 1). The received signal may be analyzed in both the 
frequency and time domains. In theory, the extent and nature of damage 
in the region of material between the sending and receiving trnsducers 
can be determined. Vary [1] has developed a counting scheme to perform 
this analysis in the time domain while more recent techniques have 
attempted to correlate various moments of frequency spectrum or the 
received signal with stiffness and strength [1,11,15]. Vary and 
subsequently other investigators have termed this quantitative parameter 
the Stress Wave Factor (SWF). A major limitation of the current 

technique is that this parameter changes dramatically with changes in 


65 



Proctor Transducer 5 Mhz Ultrasonic 

,, * Transducer 
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Procedure for calculating the SWF measure. 



such experimental variables as the receiving and sending transducer 
specifications, the couplant used, gating, sample interval and threshold 
voltage. Thus, the SWF, as it exists now, yields only information on 
the relative strength, stiffness, and damage between different specimens 
tested with the same parameters. A current research interest is to 
improve this scheme so that it will yield a meaningful quantitative 
value independent of the previously mentioned conditions. Separating 
the signal modifications caused by the material from those resulting 
from the interfaces and the electronics is an important step in this 
process [4,8], Fig. 2. After this material transfer function is 
isolated, a greater understanding of wave propagation in composites will 
be made possible and the transfer function may then be related to the 
precise state of damage in the material. 

Of primary importance in understanding the manner in which 
ultrasonic waves interact with the material is the determination of the 
type of wave that is propagating. An extensive literature review 
yielded no information on the nature of wave propagation in composite 
plates. Literature on wave propagation of anisotropic materials focuses 
on such factors as grain size and orientation which have little direct 
relevance to propagation in fiber-reinforced composite material. 
Fortunately, there is an abundance of literature on wave propagation in 
isotropic materials, some of which is relevant to composites. 

The two predominant wave types involved with ultrasonic analysis 
are Rayleigh and Lamb waves. These waves vary in the way in which they 
propagate and in the way they intereact with the material [6,7 ], Fig. 3. 

Rayleigh waves are elastic waves which propagate near the free 
surface of a solid and decay rapidly with depth. Because Rayleigh waves 
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I 

MATERIAL TRANSFER FUNCTION 
I 

INTERFACE TRANSFER FUNCTION 
1 

RECEIVING TRANSDUCER TRANSFER FUNCTION 

SIGNAL CONDITIONING TRANSFER FUNCTION 

1 

Output 


Fig. 2 Progressive altering of input prior to final output. 
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Direction of Wave Travel 



Direction of Wave Travel 



Direction of Wave Travel 



Fig. 3 a.) Symmetric Lamb wave motion, b.) Assymmetric Lamb 

wave motion, c.) Surface wave motion 
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spread in only two dimensions, they fall off more slowly with distance 
than do other three-dimensional wave types. This geometric decay is 
much less than that caused by scattering. Rayleigh waves interact very 
strongly with surface defects which causes significant attenuation. 
Absorption is analogous to that seen in longitudinal and transverse 
wave's because the Rayleigh wave is a combination of these two wave 
types. Rayleigh waves do not exhibit phase velocity dispersion. 

Lamb (or plate) waves propagate in a plate with free boundaries. 
Displacements occur both in the direction of wave propagation and 
perpendicular to the plane of the plate. Lamb waves exhibit dispersion 
and tend to interact more uniformly through the volume of the 
material. These waves tend to occur if the wavelength is greater than 
or equal to the thickness of the plate. They have, in general, both 
symmetric and antisymmetric components. These two modes propagate 
independently of one another. 

The purpose of this investigation is to use an improved transducer 
developed at the National Bureau of Standards in the study of wave 
propagation in composite materials and thus improve the state of the art 
in acousto-ultrasonic testing. It is also hoped that the use of this 
transducer will yield some insight as to the wave type being propagated 
in this technique. 

In the past, most AU work has been done using conventional 
piezoelectric pulse/echo transducers. These transducers have been 
designed for both reception and transmission of ultrasonic signals. In 
pulse/echo techniques the signal is travelling parallel to the 
transducer axis and the transducer is designed with this in mind. 
However, propagation in All techniques is generally perpendicular to the 
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axis, presenting far different demands on the transducer. Most of the 
pulse/echo transducers have a narrow bandwidth, making them relatively 
insensitive to part of the wide spectrum of frequencies in AU 
analysis. In order to further develop AU techniques, these limitations 
must be overcome through the use of an improved transducer. 

For AU applications, a transducer must have a flat frequency 
response over a wide frequency range. It must produce electrical 
signals that are accurate depictions of one independent physical 
quantity, such as normal surface displacement, while showing negligible 
response to other parameters such as velocity, acceleration, and 
tangential displacement [8]. Finally, the ideal transducer would have a 
very small aperture, thus limiting the effect of surface integration and 
the interference effects associated with signal direction. One such 
transducer recently developed by Proctor at the National Bureau of 
Standards is the NBS Conical Transducer. 

The NBS Transducer [9] is a piezoelectric transducer consisting of 
only two pieces, a truncated cone-shaped piece of lead-zirconate- 
titanate, which is the active element, and a brass cylinder backing. 
The 1.5 mm diameter of the truncated cone surface is quite small in 
relation to wavelengths of the signals of interest and in relation to 
the 13 mm diameter transducers most often used previously. This small 
contact area size makes the transducer a close approximation to a point 
receiver. The backing of the transducer serves as the non-grounded 
electrode for signal detection. Figure 4 illustrates the geometry and 
structure of the transducer. 

Glass capillary break tests have documented the NBS transducer's 
ability to electronically depict surface displacements accurately. A 
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Schematic diagram of the NBS conical transducer. 

(after Proctor) 


step function input is applied to a large thick steel specimen by 
breaking a glass capillary on the surface. Fig. 5. In Figs. 6 and 7, 
the electrical signal received a short distance away by the NBS 
Transducer is compared to the theoretical surface displacement. As 
shown, the NBS Transducer signal is almost exactly as predicted by 
theory [10]. 


II. Experimental Procedure 

The design of the conical transducer is such that the tip of the 
transducer must be grounded. This creates some problem with use on 
composite specimens due to the relatively poor conductive epoxy resin on 
which the transducer is to rest. This was thought to be remedied 
through the application of a thin coat of conductive epoxy paint (Bipox 
Traduct). This paint was applied to the region of the 0.5" wide, 
6" long graphite epoxy test specimen to be contacted by the transducer 
and an electrical lead was attached. The voltage drop across the region 
was checked and found to be minimal. The transducer was then placed on 
the specimen. The ground wire was attached to the shielding which was 
in turn wired to an aluminum base plate upon which the specimen rested. 

Electrical noise was the dominant feature of the output signal. It 
was thought that possibly the reason for the poor response of the 
transducer was lack of out of plane surface displacements in the wave 
type being propagated in the material. It has been reported by Henneke 
and Stiffler that acoustic emission signals in composites are primarily 
in-plane in nature. Thus, if AU signals propagate in the same manner as 
acoustic emission signals generated inside the specimen, poor response 
would be expected. Support for this argument was generated when a PVDF 
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Fig b Schematic of the surface-pulse and through-pulse 

(epicenter) experimental arrangements, (after Proctor) 





Fig 7 The seismic surface pulse according to Pekeris 

Abscissa, time in microseconds. Ordinate, normal surface 
displacement in picometers. Input, point-force-step-function 

(after Proctor) 



transducer (sensitive primarily to in-plane displacements) was attached 
to the specimen and a strong signal was generated. However, when the 
sending transducer was placed on the underside of the specimen and a 
situation was created in which out of plane displacements were without a 
doubt present, no signal could be received by the NBS transducer. Other 

reasons for this problem were then sought out. 

Electrical noise continued to be the only feature of the output 
when inputs were from lead break tests, from glass capillary break 
tests, and from transducers ranging from 5 to 15 mHz. Couplants with a 
wide range of viscosity were also used. Additional amplification was 
attempted as were alternate triggering schemes. The transducer was then 
rested on an aluminum plate and a strong clear signal was received with 
all of the above mentioned inputs. Thus, although at first it appeared 
to be an adequate method of grounding, the unpredictable electrical 
properties of the graphite epoxy composite prohibited it from being 

grounded with the conductive paint. 

Because the grounding on the aluminum plate seemed to be adequate, 
the conductive epoxy was abandoned in favor of a thin piece of aluminum 
foil. The foil was bonded to the composite specimen with a thin layer 
of Panametrics couplant. The other end of the aluminum foil was 
attached to the transducer shield. It was hoped that the foil would be 
thin enough so as to not affect the signal, yet thick enough to allow 
adequate grounding. A strong and relatively noise free signal resulted 

from all of the signal input methods. 

It was decided that in order to minimize reflections and edge 
effects and better isolate the AU signal of interest, a wide graphite 
epoxy plate should be used for the remainder of the investigation. A 
unidirectional 25-ply, .125" thick, 12" square plate was used. A grid 
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IMPLICIT REAL *8 (A-H,0-Z) 
INTEGER N,IWK( 14), IA( 16384) 
REAL* 8 AI(8193) , WK( 1 ) , AMP ( 8193 ) 
REAL *8 PHASE (8193 ) 

REAL *8 TIME ( 8193 ) , FREQ( 8193 ) 

900 FORMAT( A16 ) 

COMPLEX* 16 X( 8193 ) , CDABS 
READ ( 10 , 900 ) GARBAGE 
DATA PI, N/3 . 141592654,2048/ 
SMPINT=0. 1 
TTOT=SMPI NT *16384. 

NC0UNT=0 

READ(10, 1000) (IA( I) , 1=1, N) 

DO 5 MM=2049, 16384 
5 I A ( MM ) =0 
AIMAX=0 . 

TMXAMP=0 . 

N=16384 


DO 10 K=1 , N 
AI (K)=IA(K) 

IF(AI(K) .GT.127. )AI (K)=AI (K)-256. 
IF(DABS(AI (K) ) . GT.AIMAX)AIMAX=DABS(AI (K) ) 
10 CONTINUE 


20 

30 


1000 

1010 


CALL FFTRC ( AI , N, X, IWK,WK) 

NNH=N/2 + 1 
DO 20 K=1 , NNH 
TIME(K)=K*SMPINT 
AMP ( K ) =CDABS ( X ( K ) ) 

FREQ ( K ) =2 . * P I * K/TTOT 
TN=AIMAG(X(K) )/REAL(X(K) ) 

PHASE ( K ) =DATAN( TN )* 180 . /PI 

IF ( PHASE (K) . LT.O. ) NCOUNT=NCOUNT+ 1 

PHASE ( K ) =PHASE ( K ) +NC0UNT*90 . 

I F ( AMP ( K ) . GT . TMXAMP ) TMXAMP=AMP ( K ) 
CONTINUE 
DO 30 K= 1 , NNH 
A I (K)=AI (K)/AIMAX 
AMP ( K ) =AMP ( K ) /TMXAMP 
WRITE (20, 1010) (AI ( I ) ,TIME(I) , 1=1, 1024) 

WRITE ( 30, 1010 ) (PHASE (I), FREQ( I), 1=1,350) 
WRITE (40, 1010) ( AMP ( I ),FREQ(I), 1=1,200) 
FORMAT ( 13,1714) 

FORMAT (2E13. 5) 

END 


Fig. 8 Fortran program used to compute the FFT. 
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Of radial lines eminating fro™ the transducer at 10 deg Increments was 
drawn on the plate. Quarter circles of 2" and 4" radii were then drawn 
so that a matrix of 20 data locations was generated. Fig. 9. Thus the 
effects of distance from the source and angle relative to the fiber 

direction could be analyzed. 

The AU signal was generated by a Panametrics pulser and sent into 
the plate via a 5 MHz Panametrics wideband acoustic emission 
transducer. The transducer was coupled to the specimen with a high 
viscosity Panametrics couplant. Following reception by the Proctor 
Transducer, the signal is sent through a preamplifier. The signal is 
then digitized, stored, and sent to an oscilloscope for viewing. The 
Biomation signal digitizer was triggered by a signal from the pulser. 
The settings on the Biomation were as follows: 

Input Range 1 volt 

Delay Following Trigger 2-4 y sec 

Sample Interval °- 05 u sec 

Trigger Level 0,5 volts 

When the distance between transducers was increased to 4”, the 
delay was increased to 2.9 „ sec because of the longer distance that the 
signal had to travel through the material. The delay also had to be 
increased (to 3.1 u sec) after the 50 deg reading corresponding to the 
slower velocity associated with propagation perpendicular to the 
fibers. Subsequent tests were run with sample intervals of 
0.1 v sec and 0.2 u sec and without preamplification. It was determined 
that the preamplification had no effect on the nature of the signal and 
that the original instrumentation readings generated the best results. 
Therefore, all analysis was performed on the original 20 data sets. 
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Fig. 9 Schematic diagram of transducer locations 
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The digitized signal was processed by the Virginia Tech mainframe 
computer system. The program, shown here, generates plots of amplitude 
versus time, and using a Fast Fourier Transform algorithm, generates 
plots of phase angle versus frequency and amplitude versus frequency. 
These plots were then analyzed to determine the behavior of the wave in 

the material . 


HI. Discussion and Results 

Most of the current research in acousto-ultrasonics is concerned 
with analysis of the frequency spectrum of received signals. It is not 
then surprising that the most interesting results of this investigation 
are the plots of signal amplitude versus frequency. However, there are 
a few significant points to note concerning the phase-frequency and 

amplitude-time graphs. 

The primary importance of the amplitude-time plot is that it 
consists of the raw data from which all other information is 
generated. The gating, threshold, and other instrumentation settings, 
which have a great effect on the subsequent signal analysis, are 
determined by looking at the amplitude-time curve on the oscilloscope 
and adjusting the instrumentation to improve the signal quality. For 
example, in all of the curves we see a spike of maximum amplitude within 
10 to 15 microseconds, followed by a decay of the signal. Because the 
Biomation is triggered by the pulser and we wish to keep most of the 
decaying signal within our range of data collection (2048 points * 
0.05 y sec/point=102.4ey sec) the delay time had to be increased 
corresponding to slower velocity perpendicular to the fibers and longer 
distance form the 4" points. Because the plots are normalized and the 
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delay time was changed between the 2" and 4" data sets, the hard-copy 
plots shown here do not show that the signal amplitude decreased as the 
distance from the NBS Transducer was increased. Neither do these plots 
show that the maximum amplitude spike occurred later in the cross-fiber 
direction, although necessary adjustments in instrumentation indicate 
this. The plots do indicate broadening of the curves as the angle 
increases with respect to the fiber direction. Figs. 10,11. Thus, in 
the direction of least stiffness we have a flat, broad signal in 
contrast with the sharp burst associated with the fiber direction. This 
broadening indicates that dispersion is likely to be occurring. 

Dispersion occurs when different frequencies propagate at different 
rates. The incident pulse is a sharp burst and Fourier analysis 
indicates that sharp bursts contain many frequency components [12]. In 
a medium that exhibits dispersion, the amplitude versus time curve 
broadens due to these differences in velocity between the frequency 
components. 

The dispersion of the signal is indicative that Lamb type waves may 
be propagating. Other time domain evidence of the existence of plate 

wave modes is the fact that both in-plane and out-of-plane displacements 
were recorded [18]. 

The phase angle versus frequency plot (Figs. 12a, 12b) shows a very 
consistent linear relationship between phase angle and frequency. 
Previous investigators [13] have shown that this linear relationship is 
characteristic of unflawed materials, correlating sharp changes in 
elastic modulus with changes in the derivative of this plot. The use of 
phase angle data is still in its beginning stages. This recent 
development is a direct result of improved transducers, such as the 
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Fig. 12b. 


Phase Angle vs. Frequency for 0 and 90 Degrees 








conical transducer, because these transducers enable one to separate 
material prompted slope variations from those created by the integration 
process over the finite dimensions of the transducer surface. 

It also must be noted that in all cases, the slope of the 4 data 
plot was 5 to 10% higher than that from the shorter distance at the same 
angle. No trends were noticed in comparing slopes from various 
angles. Further analysis and development of the technique must be 
conducted before any conclusive statements regarding the significance of 
this phase angle information can be made. 

There are four major characteristics of the plots of amplitude 
versus frequency graphs. Figs. 13-22. They are: 

— spreading of the frequency spectrum as the angle approaches 90 
degrees, high amplitudes are more likely to exist at higher 
frequencies in the cross-ply direction. 

regular spacing between peaks and valleys of all frequency 

spectra 

increased number of peaks and valleys within a given frequency 

band as the angle approaches 90 degrees, more narrow peak and 
valley width at large angles 

-- received frequencies were predominantly in the range of 0.2 to 
2.0 MHz 

It is easy to trace the development of the SWF technique which is 
based on moments of the frequency spectrum by noting the spreading of 
the spectrum due to the increased stiffness. For example, in the 0 deg 
direction, no amplitudes greater than 0.15*V(max) exist at frequencies 
greater than 1.4 MHz, whereas in the 90 deg direction these amplitudes 
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Fig. 16. Amplitude vs. Frequency for 30 Degrees. 
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exist out to 2.0 MHz. This trend is continuous through the angles 
tested and in both 2" and 4" data sets. This spreading tends to 
indicate an inverse correlation between stiffness and first moment of 
area of the power spectrum. 

The regular spacing between the peaks and valleys of the spectrum 
can be explained by bandstop and bandpass phenomena which are 
characteristic of plate waves. Selected frequencies interfere 
destructively upon reflection from the back of the plate before they can 
propagate down the plate to the receiving transducer. These frequencies 

are given by [ 14, 151 : 
n-nV 

w " b 

where: v = velocity of the wave 

b = plate thickness 
n - some integer value 

For example, taking the thickness to be 0.33 cm and the distance 
between frequency valleys to be 0.2 MHz, the calculated velocity is 
2.1*104 cm/s. This calculated velocity is far less than the through- 
the-thickness velocities cited in the literature for graphite epoxy 
laminates of 2.0*10~§ to 18*1(T§ cm/s. This numerical inconsistancy can 
best be explained by differences in velocity and reflections from matrix 
rich layers between plies, the thickness of the couplant, and from the 
thickness of the foil layer. 

In light of previous SWF work done, it was encouraging to note that 
the predominant frequencies were below 2 MHz. A previous investigation 
[171 yielded similar results and there was some doubt as to the validity 


1-99 


of the tests due to a problem in the FFT algorithm and due to the 
transducer used. This data seems to substantiate the previous work 
despite those problems. The reason for the predominance of the low 
frequencies despite the fact that the incident pulse was centered around 
5 MHz may involve the resonant frequency of the plate although a 
different plate geometry was used 'in the previous investigation. 

It is thought that most of these characteristics can be explained 
by considering the possible conversion of the incident longitudinal wave 
to a plate type wave. This process is quite complex and much of the 
mechanism is still not clearly understood. However, it is thought that 
the signal which is incident normal to the plate begins to propagate 
along the plate due to a difference in pressure between the transducer- 
plate interface the the air-plate interface or through the excitation of 
resonance frequencies of the plate. The plate, in addition to its 
incident z component of velocity, begins to develop both x and y 
components. The magnitude of the x or y component compared to that of 
the z component, and thus the number of rejections from the front and 
back surfaces necessary to travel a given distance, is dependent on the 
orientation of the fibers relative to that direction. The recorded 
frequency spectrum is a function of both number of reflections and the 
manner in which the frequency components interact upon reflection. This 
indicates the complex manner in which the frequency spectrum depends on 
the fiber orientation [18]. Much more analysis is needed before this 
understanding is clear. 
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IV. Conclusions 


The NBS Conical Transducer has been found to be a valuable tool in 
the analysis of wave propagation in composite materials. This 
illustrates the possible role of this type of transducer in the 
improvement of the SWF scheme. Future modifications of the transducer 
should include an alternate scheme of grounding so that the transducer 
can be conveniently applied to nonconductive materials. Frequency 
spectrum analysis coupled with information from the time domain has 
substantiated the idea that plate wave modes are being propagated in the 
Stress Wave Factor analysis of thin composite plates. 
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Abstract 


The quality assurance and mechanical serviceability of advanced 
filament reinforced composites is of great concern because of the 
increasing use of these materials in many engineering structures. For 
quality control and maintainability of composites, tests based on strong 
fundamental and well understood physical principles must be established 
which can be used to ascertain and predict the design serviceability of 
finished composite parts. The acousto-ultrasonic, or stress wave 
factor, method has been found to be especially useful for providing a 
quantitative parameter which correlates well with certain mechanical 
properties of composite materials. However, our basic understanding of 
the type of wave modes established in the material when one applies the 
SWF technique is insufficient. This work investigates the propagation 

of Lamb, or plate, waves in anisotropic, homogeneous, materials. By 

( 

influence of analytic and experimental observations, it is suggested 
that the stress wave factor technique causes the establishment of the 
first few fundamental plate wave modes in a composite laminate. It is 
in large part propagation characteristics of these plate modes which an 
investigator measures when applying the SWF method. 
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I. INTRODUCTION 

The quality assurance and mechanical serviceability of advanced 
filament reinforced composites is of great concern because of the 
increasing use of these materials in many engineering structures. 
Because of the inherent anisotropy and inhomogeneity of composite 
materials, damage development and failure is distinctly different from 
that of homogeneous isotropic materials. Damage can no longer be 
defined as the formation and growth of a self-similar crack until a 
critical flaw size is reached, at which time the crack catastrophically 
propagates under the load condition. Rather, damage development in 
composite materials is a complex sequence of matrix cracks, 

delaminations, fiber breakage, fiber-matrix disbonds, and so forth, 
which form and coalesce until some critical damage state is reached, at 
which point the structure fails. 

Because composites are a relatively new material, there are not 
many years of experience upon which engineering decisions can be 
based. There is no simple accept-re ject criterion for in-service 
structures constructed from composite materials. Composite materials 
have no analog to homogeneous isotropic materials whose critical flaw 
size can be calculated from fracture mechanics and whose flaw size in a 
structural element can be used as the accept-reject criterion. 

For example, one type of in-service damage that is not readily 
detected is impact damage. If a piece of sheet metal is impacted with 
an object having enough energy to cause damage, the damage is usually 
manifested as a visible dent in the material. However, a composite 
structure can suffer significant impact damage and can have extensive 



internal damage such as delaminations and matrix cracks, yet show no 
external visible signs of damage. 

Another way in which composites differ from other materials is that 
of fabrication. Because composites are typically manufactured in one- 
step processes, there are no intermediate stages at which one is able to 
check the quality of the material used nor for manufacturing defects 
which include porosity, matrix rich areas, matrix poor areas, improper 
curing, missing plies, misoriented plies, delaminations, interlaminar 
inclusions, etc. So, both in quality control and maintainability there 
must be established tests based on strong fundamental and well 
understood physical principles which can be used to ascertain and 
predict the design serviceability of finished composite parts. 

A number of nondestructive test methods have been developed for the 
study of composite materials and have their own particular strengths and 
weaknesses. Experience with the thermoset matrix composites has 
indicated that generally a more complete understanding of the damage 
development and damage state is obtained if several NDT techniques are 
used. In particular, different NDT methods have been found to be more 
sensitive to different types of damage. Some methods such as X-ray 
radiography, thermography, vibrothermography, and ultrasonic C-scan give 
an indication of the field of damage. Of these. X-ray radiography shows 
the most microstructural detail, and thermography, the least. On the 
other hand, thermography gives an indication of the integrated effect of 
the damage field. Some NDT methods such as ultrasonic attenuation, 
acousto-ultrasonics, acoustic emission, stiffness reduction measurement, 
and thermography give real-time information. Some are directly related 
to the mechanical properties such as monitoring stiffness reduction and 
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measuring ultrasonic wave speed. Hence, several techniques, each 
providing complementary information, are usually required to obtain an 
understanding of the response and damage development in composite 
materials. 

Although there are a number of tests to interrogate a composite 
material; at this time there is still a need for simple tests which can 
yield quantitative parameters descriptive of the mechanical state of the 
material. In particular, one would like to have test methods which 
yield parameters which can be used to predict the important mechanical 
properties of strength, stiffness, and life. The measurement of 
ultrasonic wave speed is directly related to material stiffness and 
hence is a quantitative parameter which yields information directly 
about one of the needed mechanical properties. It may also be related 
indirectly to information concerning strength and life, although 
knowledge concerning the reasons these materials fail is in too immature 
a state at present to say anything definitive about this. 

Much work has been performed using through-the-thickness ultrasonic 
attenuation and velocity measurements to monitor damage and evaluate 
material properties in composites. However, little or no information 
can be obtained about the in-plane properties. Vary's stress wave 
factor or acousto-ultrasonics [1-3] is one technique that does give in- 
plane information. The stress wave number has been correlated with a 
number of material parameters but there presently is little 
understanding of why it works. 

It is well known that an acoustic or electromagnetic waveguide will 
allow only propagation of a wave having a particular frequency or any of 
its harmonics for a given wave vector. For the case of an acoustic wave 
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guide, the free surfaces cannot support a normal stress; thus acoustic 
waves propagating through a wave guide that have normal stress 
amplitudes other than zero at the free surface will be rapidly 
attenuated due to destructive interference. However, there will be one 
fundamental wave length or an integral multiple of this wave length that 
can propagate through the wave guide for a given wave vector and still 
maintain the stress free boundary conditions. 

Wave propagation in plates can be separated into several regimes 
depending on the wavelength with respect to the thickness. If the 
wavelength is much larger than the thickness of the plate, then low 
frequency plate modes are set up. For the low frequency extensional 
modes, the wave vector is in the plane of the plate and there is no 
variation in displacement through the thickness. This mode can 
propagate for all frequencies, is nondispersive, and can be easily 
applied to angle-ply laminates. Its interpretation is fairly 
straightforward. This may be the mode that can be most easily applied 
to the nondestructive testing of plates. 

Recently, analytical and experimental research has been conducted, 
under the sponsorship of General Electric Co. [1], on low frequency wave 
propagation in anisotropic plates with application to material 
characterization. It included: (1) a theoretical study conducted to 
compute low frequency plate phase velocities, and (2) measurement of 
the phase velocities in the principal directions of the plate. The 
variation of phase velocity with frequency was used to construct the 
fundamental, low frequency dispersion curves. 

The results of this study provide significant new insight into the 
nature of the wave propagation that occurs in the stress wave factor 
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technique. This wort has included calculating the low frequency 
fundamental plate wave modes in composite materials. It should be 
pointed out that for the analysis described here, the composite plate 
has been treated as a macroscopical ly homogeneous, but anisotropic, 
material having the elastic constants of the composite material in 
question. This assumption is reasonable as long as the wave lengths are 

in comparison with the fiber diameters, fiber spacing, and ply 
thicknesses, which is a reasonable assumption for frequencies less than 
those normally used for ultrasonic investigation (i.e., less than 10 

MHz). 

Preliminary experimental work has shown that the flexural phase 
velocity in the plate varies as a function of frequency due to the 
highly dispersive nature of the plate. The dispersion curves, i.e. 
curves for frequency versus velocity, are characteristic for a plate 
having a particular thickness and elastic properties. Once established, 
this analysis ailowed good correlation of experimental dispersion 
relations with plate extensional stiffnesses, A u and A 22 and plate 
bending stiffnesses 0 U and 0 22 . Analytical and experimental results, 
to date, show good correlation at the longer wave lengths. It appears 
from these results that fundamental plate wave modes are established by 
placing a transducer flat upon the surface of a plate. These plate wave 
modes propagate along the plane of the plate and can be detected by a 
similarly placed transducer "downstream" of the transmitter. 
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extensional and flexural wave propagation 

IN COMPOSITE PLATFS 

In order to make clear the implications of our previous analytical 
work [1], a condensed presentation of the analysis is briefly reviewed 
in this section. The low frequency approximation described earlier 

gives several simplifying factors for formulating wave propagation in 
plates, namely: 

1) For extensional and in-plane shear waves, the displacements 
are uniform across the thickness of the plate and out-of-plane 
inertial effects are negligible. 

2) For flexural waves, rotary inertia and transverse shear are 
negligible. 

With these approximations, classical laminate theory can be used as the 
basis for developing the low frequency wave equations. 

Again, it should be noted that this two-dimensional development is 
valid only for wavelengths larger than the thickness of the plate. For 
short wavelengths, the thin plate assumption is no longer valid because 
the displacements will be a function of the thickness coordinate 
necessitating an analysis using three-dimensional elasticity. 

2.1 In-plane Wave Propagation 

For in-plane waves, the stress free surface 
immediately satisfied by assuming a plane state of 
equations of motion for a plane state of stress are: 

30 1 + a °6 - 

3Xj_ sx 2 = pu l (2.1a) 


condition is 
stress. The 
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(2.1b) 


80g 

3X ^ 



pu 2 


For a laminate that is symmetric with respect to its midplane, the 
stress-strain relationship is given by 


°1 A 16 A 12 

ct 2 = 3 A 12 A 22 
°6 A 16 A 26 

where A i -'s are elements of the extensional stiffness matrix. The 
calculation of the extensional stiffness matrix is discussed in a number 
of books [4] on composite materials and is not discussed here. 

The strains are given by 


A 16 e l 

A 26 e 2 (2 ' 2) 

A 66 y 6 


e l " axj_ 

3U 2 

£ 2 ~ 

3U ^ 3U 2 

Y 6 ~ 3X 2 + 3*^ 

The displacements for a plane harmonic wave 
of the plate are: 


travelling 


( 2 . 3 ) 


in the plane 


u, = Aa^ exp (i (k(x 1 s. 1 +x 2 «, 2 )-u>t) 
u 2 = Aa 2 exp ( i ( k ( x ^ 5, ^+x 2 ?, 2 ) -uit ) 


( 2 . 4 ) 
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where 



A = amplitude of displacement 

“i = displacement component in the x^* coordinate direction 
i = /~A 

k = wave number = 2tt/\ 

= direction cosine of the wave vector in the x- 
coordinate direction 
w = circular frequency 
t = time. 

When Eqns. (1) to (4) are substituted into one another and common 
terms cancelled the following set of homogeneous equations is obtained. 


/^Ai 6+ *i*/i6+^A 6 6-p(£) 2 

\ l l A 16 +4 l t 2< A 12 +A 6 6 ) +l 2 A 26 


h A 16 +a lM A 12 +/ W + 4 A 26\ l a l 

a l A 66 +2s 'l S '2 A 26 +a 2 A 22' p ^^7 \ a 2; 


= 0 (2.5) 


For a nontrivial solution, the determinant of the 2x2 matrix above 
must be equal to zero. This gives a characteristic equation whose 
eigenvalues w /k are the phase velocities in the plate. The determinant 
has two roots which correspond to a quasi-extensional and quasi-shear 
wave. The respective eigenvectors and a ? can then be obtained from 

Eqn. (2.5). The wave velocities can be determined for any direction in 
the plate. 

For the experimental work, the velocities were measured in the x^ 
direction and the direction and only these directions will be 
considered in the subsequent analysis. 

For 8^-1 and 8 ,, = 0, the extensional wave velocity is given by 
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_ 4 

V 1 = ^ (A U +A 26 ) + / (A n +A 66 ) 2 + 4[A n A 66 -A 16 2 )^ 

and 


a l 

a 2 


= 1 
= 0 


That is, for this propagation direction the wave is purely longitudinal. 

For a unidirectional laminate or a laminate symmetric with respect 
to the x-l axis, the extensional velocity is given by 


Likewise for propagation along the x 2 axis, a- 2 1 an< ^ ^ 0 gives 

_ h 

V 2 = (0 < A 22 +A 66> + / (A 22 -A 66 ) 2 + KAjjAge-Aje)' <2 ' 8) 


and 

1 
2 


For a unidirectional laminate or a laminate symmetric with respect 
to the x 2 axis the in-plane shear wave velocity is 


M 


2 



(2.9) 


2.2 Flexural Wave Propagation 

The derivation for flexural wave propagation is carried out in a 
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similar fashion as extensional wave propagation. The governing equation 
for flexural wave propagation is 


3 2 M. 29 2 M c 

-4 + + 


3 2 M. 


9X 


1 


9X^3X 2 3X 2 


d 3?u 3 

D n 2 
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where 


M-j = bending moment about x i coordinate axis 
p = mass density 
d = thickness of plate 

The bending moments are related to the curvature by 


M 1 

D 11 

°12 

°16 

k i 


M 2 . 

°12 

°22 

D 26 

k 2 

(2.11) 

M 6 

°16 

°26 

°66 

k 6 


where D 1 j are the bending 

stiffnesses and 

k .. are 

the curvatures 

of the 

plate. The curvatures are 

given by 






K l = ~ 

9X 2 





< 2 = ~ 

d 2 u 3 

9X^ 



(2.12) 



9 2 U 3 





gx^ax^ 


The displacement of a flexural wave, u^, is given by 
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i (k ( fi. , x , + -wt] 

u 3 = Ae 11 * d 


(2.13) 


where the parameters have the same definitions as for Eqn. (2.4). 

Substituting Eqns. (2.10) through (2.13) into one another gives 


kV?D u 


+4j,^2 D 16 




+2D 66) 


4z, 1 £ 2 D 26 + 


£ 2^22^ = pdw 2 (2. 14) 


For flexural waves propagating along the x^ axis, 
reduces to 

4 ? 

k D. , = pdw 

Since 

V = — 

P k 

the phase velocity is 




pd 


4 7 r 2 D 11 f2 

■ od 


k 


Eqn. (2.14) 

(2.15) 

(2.16) 


(2.17) 


It should be noted that the phase velocity is no longer constant but is 
a function of the frequency. Thus, propagation of low frequency 
flexural waves is dispersive. 

The group velocity is 


and hence 



(2.18) 


= „ JL r£ d* ^ 

g pd2u pdw ' D j. 
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2.3 Calculated Velocities for Various Materials 

In order to verify that the analysis presented in the previous 
section represents the characteristics of the waves which are generated 
and detected during application of the acousto-ultrasonic technique, 
both analyses and experiments were performed on two different materials, 
aluminum and graphite epoxy. The aluminum plate was 6061-T6 aluminum, 
0.065" thick x 12" x 12". The graphite epoxy panels were positive 
pressure molded from AS-4(GR)/PR288 and consisted of [0]g, [01 22 , 
[0,90] 2s , and [ 0 ,+45 ,90 , -045 1 s laminates. The mechanical properties of 
these materials can be found in Table 1. The laminate properties for 
the plates can be found in Table 2 and the low frequency extensional and 
flexural wave velocity can be found in Table 3. 


Table 1. Material Properties 


6061-T6 

Aluminum 

AS-4(GR)/PR288 
Graphite Epoxy 


E 11 

(MSI) 

E 22 

(MSI) 

G 12 

(MSI) 

v 12 

V 21 

P 

(lbf-sec^/in 4 ) 

10.00 

10.00 

3.75 

0.30 

0.30 

1 . 55x 10 -4 

18.50 

1.30 

0.67 

0.32 

0.02 

1. 50xl0 -4 
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Table 2. Laminate Properties 


Material 

Lay-up 

Thickness 

(inches) 

Q 11 

(MSI) 

Q 2 2 

(MSI) 

A 11 

JJlxlO 6 

a 22 

JL^xlO 6 

D 11 3 

in-lbsxl0 J 

°22 

i n- 1 bsxlt 



in 

in 

in 

in 

6061-T6 









Aluminum 

— 

0.0625 

11.00 

11.00 

0.69 

0.69 

0.224 

0.224 

Graphite 

Epoxy 

1 0 1 8 

0.092 

18.6 

1.31 

1.71 

0.12 

1.21 

0.085 

Graphite 

Epoxy 

[ 0)22 

0.262 

18.1 

1.31 

4.87 

0.34 

27.8 

1.96 

Graphite 

Epoxy 

[0,90] 2s 

0.095 

— 

— 

0.940 

0.940 

.942 

.480 

Graphite 

Epoxy 

[0,45,90,-45] 

l s 0.097 



0.747 

0.747 

.891 

.268 
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Table 3. Low Frequency Extensional and Flexural Wavelengths 


Extensional Waves Flexural Waves 


Material 

Stacking 

Sequence 

V 1 

(in/sec) 

V 2 

(in/sec) 

V 1 

(in/sec) 

V 2 

(in/sec) 

6061-T6 

Aluminum 

— 

207,695 

207,695 

153/T 

153/” 

Graphite 

Epoxy 

[0] 8 

352,136 

93,452 

242/~f 

125/” 

Graphite 

Epoxy 

[ 0)22 

352,136 

93,452 

405/T 

210/” 

Graphite 

Epoxy 

[ 0 ,90 ] 2 S 

256,836 

256,836 

226/T 

191/” 

Graphite 

Epoxy 

[0,45 ,90,-45 ] s 

226,583 

226,583 

2 23/7 

164/” 


11-17 


III. EXPERIMENTS 


A number of experiments utilizing low frequency acoustic waves were 
performed on several different laminates and aluminum to determine the 
applicability of the plate wave analysis presented in the previous 
section. A piezoelectric transducer, which acted as the acoustic 
generator, was excited using a gated sine wave. The gated sine wave 
provided an excitation signal with a very narrow frequency band while 
still retaining a definite leading edge. 

A second piezoelectric transducer was used as the receiver. The 
position of the receiver was advanced a known distance and the change in 
time for a particular phase point noted from a digitizing 
oscilloscope. Dividing the displacement by the time yielded the phase 
velocity, from which the laminate stiffnesses were easily obtained. 

3.1 Experimental Arrangement 

Figure 3.1 is a schematic of the experimental arrangement for 
determining the phase speed in the plate. The excitation signal was 
generated by a signal generator which allowed the signal to be gated. 
The frequency of the signal was adjusted by the signal generator. The 
gate signal was created by a pulse generator. The amplitude, width, and 
repetition rate of the pulses were adjustable. The frequency of the 
signal was obtained from a frequency counter. 

The signal from the signal generator, which had a peak-to-peak 
voltage of 20 volts was input into a piezoelectric transducer which 
converted the electrical signal into mechanical vibrations. The 
transducer was a Panametrics V401, l"x % ", wide band transducer having 
a center frequency of 1 MHz. The transducer was acoustically coupled to 
the plate by a water soluble couplant. 
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Schematic of experimental arrangement for obtaining phase 
velocities in plates. 




The acoustic signal propagating in the plate was detected by a 
second transducer which was a Panametrics V109, % inch diameter, wide- 

band transducer having a center frequency of 5 MHz. The signal was 
amplified by a preamplifier and an amplifier and displayed on a 
digitizing oscilloscope. Both the transmitting and receiving 
transducers were mounted flat against the specimen surface. 

The oscilloscope had movable cursors and a digital readout which 
corresponded to the position of the cursors in the time domain. This 
gave an accuracy in the time domain of ±0.5 microseconds. The digital 
oscilloscope had the capability of storing the waveform on disk and 
plotting the waveform on an x-y plotter. 

3.2 Experimental Procedure 

The experimental procedure consisted of first setting the signal 
generator to the desired frequency, usually between 10 kHz and 1 MHz. 
The cursor on the digitizing oscilloscope was then positioned on a 
particular point on the received wave, usually a peak, and the time 
noted. Next the receiving transducer was translated a fixed amount. 
The selected phase point on the wave was followed and the time noted 
again. The phase velocity is given by 

\/ = A 

p at 

where 

d = distance receiving transducer was translated 
at = change in time of point of phase 
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3.3 Data Reduction 


The wave number was obtained from 



for the frequencies at which data were obtained. Next the frequency and 
wave number were multiplied by the thickness of the plate which gave fd 
and kd. A plot was then made of fd versus kd. This is a dispersion 
curve, from which the phase velocities for any fd are obtained by taking 
the secant to the curve and the group velocities are obtained by taking 
the tangent of the curve at the value fd. 
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IV. RESULTS 


Numerous tests were conducted on plates made from different 
materials and having a variety of thicknesses to determine the 
plausibility of applying the plate wave analysis to the wave modes 
transmitted in thin plates when one applies the stress wave factor 
technique. While obtaining data, it was noted that as the frequency of 
the excitation signal was increased, there existed a frequency at which 
there was an abrupt change in the phase velocity. Further investigation 
suggested that this was a change in mode from flexural wave propagation 
(the lower frequency signal) to extensional wave propagation (the higher 
frequency signal). Also, from observing the received signal, it was 
noted that the signals for the lower frequencies were dispersive while 
the signals at the higher frequencies were not. 

Plots of kd versus fd were made with the values obtained from the 
analysis and from the experiments plotted on the same graph. Figures 
4.1 through 4.7 are graphs of the experimental data (points) and 
theoretical results (curves) for the materials tested. To emphasize how 
the experimental points were taken, we note the following: Data points 
were first taken by applying the lower frequencies to the transducers. 
As the frequency was increased, the data points marched up along the 
lower curve. At a critical frequency, the generated wave mode changed 
discontinuously so that suddenly a higher velocity mode was generated. 
This resulted in a point appearing at the lower end of the higher curve 
in Figs. 4. 1-4.7. With further increase in frequency, the data points 
marched up along the higher curve. 

There is an excellent correlation between the experimental data and 
theoretical results for the higher (longitudinal) mode , kd=0.5 for all 
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Low frequency dispersion curve 
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Low frequency dispersion curve along the axis for graphite 
epoxy plate. 
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Low frequency dispersion curve along the x 2 ax i s f or graphite 
epoxy plate. 
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Fig. 4.4 Low frequency dispersion curve along x^ axis for [0,90]2 
graphite epoxy plate. 
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'equency dispersion curve along X 2 axis for [0,90]2 S 
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Low frequency dispersion curve along x;? axis for 
[0 ,45 ,90 ,-45]s graphite epoxy plate. 


frequencies tested, and for the flexural mode for frequencies up 
to kd>0.5 the thin plate assumption may no longer be valid for the 
flexural mode and a better approximation is needed. Also, the flexural 
mode is noticeably dispersive and it is somewhat more difficult to 
measure the phase speed. From these plots, the extensional stiffnesses 
and A 22 can easily be found if 'the density of the material is 
known. 

The experimental data and analysis clearly indicates that at low 
frequencies acoustic energy is propagated by the fundamental extensional 
and flexural plate modes. This could explain the nature of the waves 
generated when one applies the stress wave factor technique. 

For example, Govada [5] used a 5 MHz, wide-band transducer to 
launch an acoustic wave in a composite specimen. However, a spectral 
analysis of the received signal showed a frequency content much lower 
than 5 MHz. For an undamaged [0,90 2 l s graphite epoxy specimen, Govada 
obtained a tri-modal frequency distribution with frequencies centered 
around 250 kHz, 750 kHz, and 1 MHz. 

With the exception of the fundamental extensional and flexural 
mode, plate modes exhibit cut-off frequencies for which waves having a 
frequency content below the cut-off frequency will be attenuated. The 
lowest cut-off frequency for a plate, which occurs for an antisymmetric 
mode, is 



For the [ 0 ,90 2 1 s coupon, d = 0.030" and V $ = 60,000 in/sec gives the 
lowest cut-off frequency at 1 MHz. Plate waves having a frequency less 
than 1 MHz will propagate as one of the fundamental modes. So the type 


11-30 



of plate modes observed in the stress wave technique are in all 
probability the fundamental extensional and flexural modes. 

The reason higher modes do not propagate requires further 
investigation. It may simply be that higher modes are more rapidly 
attenuated and hence become lost in the noise level when one is making a 
measurement. 
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Part III 


NONCONTACTING DETECTION IN ULTRASONIC NDE OF MATERIALS: 

SIMPLE OPTICAL SENSOR AND 
FIBER OPTICS INTERFEROMETRIC APPLICATION 

A. Sarrafzadeh-Khoee and J. C. Duke, Jr. 

Abstract 

An optical displacement sensor based on the two-beam interferometry 
of coherent light using two principal diffracted orders of a fine 
reflection blazed grating is constructed. The inherent capability of 
the interferometer for absolute calibration, its maximum sensitivity 
adjustment, and the simplicity of the configuration of its components 
are explained. 

Utilization of a single-mode optical fiber in the test arm of the 
interferometer to enhance the geometric flexibility of the optical 
system is demonstrated. 

The sensor's application in the noncontacting detection of 

stimulated ultrasonic stress waves for nondestructive evaluation of 
advanced composite materials is emphasized. The detected signal is 
digitized and recorded using an IBM PC with a plug-in data acquisition 
board. The combination of high sampling rate capability of the device, 
Fast Fourier Transformation of data, and optical detection has 
facilitated the spectral analysis of ultrasonic signals in order to 
achieve interpretation and quantification of the physical parameters of 
interest. 
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Introduction 


The acousto-optical techniques for the detection and analysis of 
ultrasonic stress waves have been of major interest to many 
investigators due to increased ultrasonic applications in nondestructive 
testing and materials evaluation 11,2]. Most physical effects are 
transformed into displacement. Different optical set-ups can be 
constructed and optimized for looking at a particular particle 
disturbance (out-of-plane/in-plane displacement) of the material under 
study. Ordinarily these optical transducers act as displacement 
sensors; therefore, any other related physical quantities can be derived 
from the displacement magnitudes. 

The main advantages for using optical sensors are that they are 
noncontacting so that no alteration of actual surface displacement 
occurs. They can be used in hostile environments, in particular for 
high temperature materials characterization. They are sensitive — 
e.g., displacements of subangstrom magnitudes can be measured. Because 
of the ability to focus the laser beam on a spot-sized region, point-by- 
point interrogation of the surface and complete characterization of the 
displacement field is readily accomplished. Following the suggestion 
made in Ref. [3], a unique optical arrangement is designed in a 
Michelson-type interferometer using a blazed diffraction grating as a 
beamsplitter and recombiner. 

System of Two-Beam Interference 

A plane reflection blazed grating is shown in Fig. 1. The grating 
is constructed by a number of densely packed right angle prisms grooved 


1 1 1-2 



Principal 

Diffracted 

Order 



Fig. 1. Single-order blazed reflection grating, (y = blaze angle, 
g = grating pitch). 
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on the surface of a plane glass block on which a thin layer of aluminum 
has been deposited for increased relfectivity. In contrast to ordinary 
(amplitude or phase) diffraction gratings where the undiffracted (zero) 
order has the stronger intensity, in a blazed diffraction grating the 
stronger intensity is shifted to some higher order. Consequently, the 
blazed grating has' a higher diffraction efficiency compared to its 
counterparts. For a given optical wavelength, the intensity of the 
diffracted beams, in general, is the function of the shape, angle, and 
depth of the grooved surfaces and of the grating surface reflectivity. 
The angle of the principal diffraction order is solely determined by the 
proper choice of the wavelength used and by grating spacing (pitch). As 
shown in Fig. 1, the angle of incidence can also be chosen so that the 
dominant diffracted order from the fine blaze grating bisects the angle 
between the dominant zero order and the incidence. Such an arrangement 
will provide a single-order blazed diffraction grating with two dominant 
orders of near equal intensity. 

As was mentioned previously, when the correct light wavelength is 
used, when the spacing of linearly groved surfaces is proper, and when 
the selection of the blaze angle is correct, the grating behaves such 
that an incident coherent beam of the laser emerges in two dominant 
diffracted orders, e.g., the zero and the first order. Referring to 
Fig. 2, consider the conditions where a) the incident angle is in a 
direction normal to the grating plane. In this case the undiffracted 
beam also emerges in the direction normal to the grating and the 
principal diffracted beam emerges at twice the blaze angle with respect 
to normal incident direction. The latter is the natural direction of 
reflection from grooved surfaces, b) The incident angle is twice the 
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blaze angle with respect to the grating normal and oriented toward the 
blaze arrow. In this case the first order beam bisects the angle 
between the zero order and the incident direction. Interestingly 
enough, this diffracted beam is along the direction normal to the 
grating plane. c) In arrangement (b), if a mirror is placed 
perpendicular' to the first order direction and if the action of 
arrangement (a) is added, a two-beam interferometric system is 
developed. The reference arm of the interferometer is the zero order — 
the natural reflection of the incident beam from grating plane — and 
the sensing arm is the first order diffraction. The sensing arm after 
reflection from the test surface is diffracted again by the action of 
the grating. Consequently, the resulting first order diffraction of the 
original first order is superimposed upon the original zero order 
diffraction and produces an interference pattern. The system of the 
unequal-path interference of these two beams is the basic optical 
component of the entire sensor. This two-beam unequal -path 
interferometer, as well as any of the Michelson-type interferometers, is 
only sensitive to the out-of-plane displacement of the motion of the 
surface under interrogation. 

Phase Displacement vs. Optical Intensity 

The test arm and the reference arm are the two interfering coherent 
beams of the optical system. The phase change in the test arm with 
respect to the reference arm is introduced due to the distance between 
the grating and the surface under study. This distance will change due 
to surface displacement fluctuation caused by the presence of ultrasonic 
stress waves. In that case, the phase displacement is written as 
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$ = 2K1 + 2K6 


where 

1 = optical distance between grating and the test surface 
K = 2n/x , optical wave number 
6 = acos wt, amplitude of surface displacement 
a) = itf, angular frequency of the oscillation. 

The factor of two in the expression for the phase displacement 
denotes the fact that light travels twice through the optical path 
between the grating and test surface. If both superposing beams of the 
interferometer have the same optical intensity (this is a reasonable 
assumption since both zero order and first order diffractions of the 
blazed grating have relatively close intensities), upon superposition, 
the interferometric intensity deviates from a minimum of zero to a 
maximum of four times the intensity of each beam (I), depending on the 
phase difference between the two beams. In general, the cosine squared 
expression for the intensity of the two-beam interference [4] is written 

as 

I' = 41(cos <t>/2) 

= 21 [ 1+cos(2K1+2Kacos wt)]. 

If the spacing between the grating and the test surface is such that 
introduces a phase of -n/2 radians, then 
I' = 21 [ 1+cosU/2+2Kacos wt)] 

= 21 - I". 

The interferometric term is 
I" = 21 sin(2KAC0S wt) . 

If the amplitude of the surface displacement is much smaller than the 
optical wavelength, s « l , the sine can be approximated to be equal 

to its argument. Thus, 
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I" - rlkAcos ut. This linear relationship between the instantaneous 
surface displacement, 6, and the optical intensity (time-averaged of the 
optical field) is the prime factor in the use of such a two-beam 
interferometer for detection. 

The variation in the light intensity due to surface disturbance of 
the test sample can be photoelectrical ly observed. The photodetector 
converts the optical power into an electrical current. For a 
photodetector with a quantum efficiency of a, the RF signal current 
produced by the incident optical intensity is 
1 S = oI“. 

The sensitivity of the detected signal is limited to the presence 
of the noise inherent in the detecting system. This would include: 
laser noise (optical) caused by fluctuation of laser power output, 
Johnson noise (thermal) caused by the load resistor of the 
photodetector, the electronic noise in the electronic signal processing 
components, and shot noice (photon) of the square-law detector, which is 
considered to be the most limiting factor in sensitivity. Other sources 
of noise, which will be taken into consideration in a later section, are 
low-frequency environmental noises (mechanical and thermal). 

The dc optical power, 21 incident on the photodetector, gives rise 
to the rms shot noise current. That is, 

'n ' 
where 

*dc = averaged optical intensity on photo-diode 
B = detection bandwidth 

e = electron charge. 
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The minimum detectable displacement amplitude is obtained when noise 
current and signal current are equal (S/N=l), i.e., i n = i s . Thus, in 
terms of output laser power, P, the maximum detectable displacement is 
written as 

6 . = f2eB/anP)*X/4u 
min v 

where ri = optical system efficiency. The theoretical minimum detectable 
displacement for a He-Ne laser, P = 10 mW, a = 0.4 W/A, B = 5MHz, 
and n * 80% is calculated to be about 0.01 A. 

Calibration and Large Phase Shift 

The optical interferometer described earlier has the same 
sensitivity as any of the Michelson-type interferometers. The output of 
the interferometer can be calibrated so that the absolute magnitude of 
the surface displacement can be determined. If the surface under study 
is subjected to large displacement excitations (6 > x/4) and the peak- 
to-peak signal output is noted, then the magnitude of the smaller 
displacements (6 < x/4) can be obtained from its peak-to-peak signal 
output. 

As was mentioned earlier, the interference of two equal -intensity 
beams has a magnitude which varies between two extreme opposite values 
depending on the relative phase difference (<t>) between the test and 
reference beams. For the interferometer to operate at the point of 
maximum sensitivity (i.e., smallest phase shift giving rise to largest 
change in intensity), the relative phase must be kept constant halfway 
between the two limiting phases throughout the operation 15]. 
Unfortunately, the presence of low-frequency mechanical and thermal 
fluctuations existing in the testing environment alters the proper 
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tuning. To solve this problem, a slight modification to the optical 
system of Fig. 2 may suffice -- the blazed diffraction grating can be 
driven by a low-frequency piezoelectric crystal at a large displacement 
amplitude so that the face of the grating translates at least an 

amount X/4. this causes the passage of phase difference through the 
optimum point of sensitivity at the time of high-frequency transient 

surface fluctuations. The frequency of the phase shifter should be low 
enough to allow the entire signal to be detected at the maximum 

operating condition. Through the use of this large phase shifter, the 

interferometer can also be calibrated. 

Fiber Optics Interferometric Application 

The use of optical fibers as a light guiding medium can be a great 
advantage in optical probing flexibility [6], In this regard, an 
attempt has been made to utilize single-mode-step-index fiber optic 
interferometry. There are special kinds of single-mode optical fibers 
that preserve the state of linearly polarized light, which is necessary 
to maintain in the interferometric application. These types of bi- 

refringent fibers have geometric flexibility and allow the test arm of 
the interferometer to be isolated from the bulky optical table. There 
are two classes of fiber optic sensors. In one application the fiber 
itself is used as a sensing device. In another application it is used 

to simply guide the optical wave through its medium. The latter has 

been incorporated into our two-beam unequal -path interferometric 
system. The optical system shown in Fig. 3 is the optical fiber link 
extension of the probing arm of the previously discussed 

interferometer. The test arm of the interferometer in the air medium 
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. 3. Schematic diagram of a flexible single-mode fiber opti 
interferometric system. 




has been coupled into the fiber medium using a microscope objective to 
focus the beam on the front flat surface of the stripped fiber. The 
second microscope objective provides the beam illumination and reception 
from the reflective surface under study. The use of an interferometric 
fiber optic sensor in this fashion allows the ultrasonic scanning of the 
reflective surface of a stationary object which is located away from the 
bulky optical bench. 

Spurious fringes may occur as a result of the interference of some 
reflected light from the front face of the optical fiber (i.e., at the 
focused laser beam-to-fiber interface). Since these unwanted fringes 
may superimpose upon the actual interference pattern, care should be 
taken to minimize this contrast limiting effect. For the maximum fringe 
visibility, the total optical path length traveled by the test arm of 
the interferometer must be an integral multiple of twice the laser 
cavity, as long as it does not exceed the coherence length of the laser 
source. Thus, by proper positioning of the fiber coupler with respect 
to the grating, and of the probing fiber coupler with respect to the 
surface under study, and by correct choice of fiber length, one can 
establish the condition under which the maximum contrast will emerge in 
the interference pattern but not in the unwanted fringes. 

Preliminary Results 

The schematic diagram in Fig. 4 shows a data processing procedure 
used in optical detection of ultrasonic stress waves. As a 
demonstrative exercise, a repetitive electronic short pulse was 
transmitted to drive a 1 MHz PZT transducer of the thickness-dilatation 
type which was coupled acoustically to one side of a unidirectional 
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Fig. 4. Schematic diagram of an NDT using ultrasound, optical detection 
and spectral analysis. 
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fiber-reinforced composite laminate. The interrogation point on the 
same side of the specimen was mirrorized for light reflectivity. Then 
three different transducer positions about 1 inch from the point were 
chosen such that the transducer could propagate plane waves parallel to, 
off-axis to, and perpendicular to the fiber axis. The resultant surface 
displacement excitations were detected by the linear optical sensor 
(arrangement in Fig. 3). After signal amplification and suitable 
filtering, the information pertaining to the amplitude and frequency of 
acoustic disturbances was obtained. Both the time-domain and frequency- 
domain amplitudes of the detected stress waves for the aforementioned 
transducer positions are shown in Fig. 5(A)-(D), respectively. The 
spectrum of the signal is not shown for the curves (C) and (D) because 
no change in spectrum has been noticed for all three of these cases. 
Note that the magnitude of surface displacement is larger in curve (C) 
than in curve (A) and is drastically reduced in curve (D). This is 
contrary to the results found in homogeneous isotropic materials in 
which the magnitude of signals for all three transducer locations would 
have been equal. 

NDT Applications 

Low-power ultrasound is one of the NDT methods which has long been 
used for materials' mechanical characterization through the 
determination of acoustic wave parameters of interest (e.g., velocity, 
frequency, amplitude, attenuation, etc.). The sources for such an 
ultrasonic regime are generated either by the materials themselves -- 
acoustic emission -- or by external means — acoustic stimulation. To 
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Fin 5 Plots of optically detected ultrasonic signals, a) Direction 
F 9 * of propagation parallel to fiber direction, b) Frequency 
spectrum of the signal shown in a). 
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quantitatively characterize materials using ultrasonic methods, a well 
characterized quantitative detection method in conjunction with a well 
characterized source of ultrasonic excitation must be used. The 
applications of noncontacting optical sensors, along with optical 
generation of ultrasound , are conducive to a more refined analysis in 
NDE. 


Conclusion 

The analysis of this optical system shows that the stress wave 
causes a change in interferometric light intensity due to surface 
disturbance. The deviation in light intensity is detected 
photoelectrical ly. The output signal is found to be proportional to the 
instantaneous surface displacements. Through suitable calibration of 
the system, the absolute measurement of surface displacements, is 
obtained. The ultimate sensitivity of these interferometric optical 
transducers is limited in part, due to the photodetector shot-noise 
current, and in part due to optical path length differences. The small 
optical wavelength gives rise to much smaller optical path length 
differences which may be less than the atomic dimension. 

The basic optical arrangement described in this paper is a simple 
one in which the reflection blazed diffraction grating is the primary 
component of the interferometer. Through a simple rearrangement of the 
previously discussed optical configuration, e.g., by reversing the blaze 
arrow orientation and selecting the zero order as the sensing arm, a 
similar two-beam interferometer is formed. In this case, the mechanical 
movement of the grating does not introduce any optical path difference 
between the first and the zero order beams of the interferometer. 
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Therefore, unlike the Michel son interferometer, this arrangement is not 
sensitive to the low-frequency noise associated with the mechanical 
vibration of the reference mirror. 

In the ordinary Michelson interferometer, only half of the laser 
power is utilized and the other half is fed back into the laser. The 
optical intensity returning to the laser might alter the laser power 
output which in turn might cause an increase in the signal noise level 
and reduce the sensitivity. The feedback problem is usually solved by 
using optical or acousto-optical isolators [7] which effectively prevent 
any light from entering the laser. In another treatment of the problem, 
a dual-quadrature Michelson interferometer [5,8] is constructed 
basically to use advantageously the complementary optical intensity 
pattern which was fed back to the laser. For our interferometer, which 
utilizes a single-order blazed diffraction grating, most of the energy 
of the incident laser beam is distributed into the zero and the first 
order so that very little light corresponding to the second order 
diffraction is returned to the laser. 

The addition of an optical fiber to the test arm of the 
interferometer as a light-guiding medium has been shown to be a great 
advantage in optical probing flexibility. For this reason, an attempt 
to utilize single-mode fiber optic interferometry for remote sensing is 
justified. 
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Closure 


Ultrasonic methods are unequaled in directly and nondestructi vely 
evaluating the mechanical behavior of materials containing damage. The 
work supported by this grant has expanded the basic understanding of the 
acousto-ultrasonic method for that purpose. 

A summary of the more significant conclusions is listed here. 

Part I 

1. Acousto-ultrasonic examination of fatigue damage development 
reveals a sensitivity to damage mechanisms that are 
undetectable by penetrant-enhanced stereo X-ray radiography or 
stiffness monitoring procedures. 

2. "Weak" regions detected by the acousto-ultrasonic method 
coincide with the location of initial damage development and 
final failure for uniaxial loaded specimens. 

3. Incorporation of the NBS (Proctor) conical transducer into the 
acousto-ultrasonic procedures has revealed that a significant 
portion of the ultrasonic energy is transferred throughout the 
laminate by plate wave propagation. 

Part II 

4. Careful analytical consideration of in-plane wave propagation 
has shown that the dominant frequencies observed in the 
acousto-ultrasonic method correspond to the fundamental 
extensional and flexural plate modes. 

5. Measurement of plate wave velocities allows for determination 
of extensional and flexural stiffnesses from place to place in 
the laminated plates. 
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Part III 


6. Practical, quantitative, noncontact detection of ultrasonic 
waves used for effecting the acousto-ultrasonic method is 
possible if optical fibers are integrated into a blazed grating 
Michel son-type interferometric system. 



1. Report No. 

NASA CR-3976 


4. Title and Subtitle 


2. Government Accession No. 


Ultrasonic Stress Wave Characterization of 
Composite Materials 


7. Authors) 

John C. Duke, Jr., Edmund G. Henneke II, 
and Wayne W. Stlnchcomb 



9. Performing Organization Name and Address 

Virginia Polytechnic Institute and State University 
Engineering Science and Mechanics Department 
Blacksburg, Virginia 24061 - 4899 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 


3. Recipient's Catalog No. 


5. Report Date 

May 1986 


6. Performing Organization Code 


8. Performing Organization Report No. 

None 


10. Work Unit No. 


1 1 . Contract or Grant No. 

NAG3-323 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 

506- 53-1 A (E- 2948) 


15. Supplementary Notes ~ “ ” “ “ 

Final report. Project Manager, Alex Vary, Structures Division, NASA Lewis 
Research Center, Cleveland, Ohio 44135. 


16. Abstract " ” 

The work reported covers three simultaneous projects. The first project was con- 
cerned with (1) establishing the sensitivity of the acousto-ultrasonlc method for 
evaluating subtle forms of damage development In cyclically loaded composite 
materials, (2) establishing the ability of the acousto-ul trasonlc method for 
detecting Initial material Imperfections that lead to localized damage growth 
and final specimen failure, and (3) characteristics of the NBS/Proctor sensor/ 
receiver for acousto-ultrasonlc evaluation of laminated composite materials. The 
second project was concerned with examining the nature of the wave propagation 
that occurs during acousto-ultrasonlc evaluation of composite laminates and dem- 
onstrating the role of Lamb or plate wave modes and their utilization for char- 
acterizing composite laminates. The third project was concerned with the 
replacement of contact-type receiving piezotransducers with noncontacting laser- 
optical sensors for acousto-ultrasonlc signal acquisition. 


17. Key Words (Suggested by Author(s)) 

Nondestructive testing/evaluation; Stress 
waves; Acousto- ul trasonlcs ; Fiber- 
reinforced composites; Damage; Strength; 
Stiffness; Signal analysis; Laser/ 
piezotransducers 

18. Distribution Statement 

Unclassified - unlimited 
STAR Category 38 

oecurny ^lassir. (or this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of pages 

158 

22. Price* 

A08 

* 

1 


For sale by the National Technical Information Service. Springfield, Virginia 22161 


NASA-Langley, 1986 





